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I. INTRODUCTION 



A. FLOW SEPARATION 

Real fluids flowing over bluff bodies separate under the influence of 
adverse pressure gradient as momentum is consumed by wall shear. Separa- 
tion is best understood in terms of its consequences. In fact, it was the 
keen observations of separation in diffusers that led Prandtl to introduce 
the concept of transition or boundary layer theory. Even though the 
boundary layer theory has revolutionized fluid dynamics, the phenomenon 
which gave impetus to its discovery still remains an enigma. 

Through the years various characteri sties of separated flow about 
bluff bodies in general and about circular cylinders in particular have 
been partially understood through observations, measurements, and analyses. 

The separation point may be mobile as on a circular cylinder or fixed 
along a line as on a wedge or blunt-based body. When a separation point 
is mobile, the pressure decreases up to the point of maximum velocity 
(slightly upstream of separation) and then increases up to the point of 
separation. When the separation point is fixed, the pressure continues 
to decrease up to the point of separation. Consequently, there are signi- 
ficant differences in the separation of flow about a circular cylinder 
and a sharp-edged body. In either case the velocity distribution, rate 
of vorticity flux, and the amplitude and frequency of oscillation of the 
separation points are not independent of the wake in the downstream side 
of the body. It is this interaction that makes the calculation of sepa- 
rated flows exceedingly complex. In other words, one needs the evolution 
of the boundary layer on the forebody (upstream of separation) of the 
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cylinder in order 1o delermine tlie i nstunidneoub position of the separa- 
tion points. But, the characteristics of the boundary layer are, in turn, 
strongly affected by the fluid motion in the wake. Evidently, one needs 
some kind of an interactive mechanism whereby the effect of the wake on 
the boundary layer, and hence separation, can be given proper considera- 
tion. In mathematical terms one can state that the unseparated potential 
flow, or the d’Alembert flow, does not even constitute a first approxima- 
tion to the separated real flow. Partly the realization of this fact, 
partly the difficulties associated with the integration of the Navier- 
Stokes equations across a singularity (separation point), and partly lack 
of sufficient computational capability have prevented the prediction of 
the kinematic and dynamic character i st ics of separated flows about bluff 
bodies. In the following, the basic characterist ics of flow about circu- 
lar cylinders are described with special emphasis on those aspects which 
prepare the reader for the understanding of the numerical model developed 
during the course of this work and the assessment of the similarities and 
dissimilarities between the predictions and observations and measurements. 

B. VORTEX SHEDDING 

The separation of flow gives rise to alternate vortex shedding. This 
phenomenon is intrinsic to the flow itself and has been known at least 
since the times of Leonardo de Vinci. In 1878, almost exactly 100 years 
ago, Strouhal discovered that there is a relationship between the frequency 
of vortex shedding, the velocity of flow, and the diameter of the cylinder. 
This relationship, which is now known as the Strouhal Number, St, is 
given by St = fD/U. Strouhal’s experiments gave an average value of St = 
0.185. During the course of numerous experiments since then, it has been 
experimentally demonstrated that: (i) the Strouhal number is dependent on 
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the Reynolds tiurnber and, as first suggested by Lord Rayleign, must be 
written as f = St(Re)U/D; (ii) the vortex shedding frequency is only an 
average quantity for a given flow and that there are always secondary and 
tertiary frequencies present; and finally, (iii) there are regions of flow 
such as the region of drag crisis where it is quife difficult to identify 
a particular frequency. It is evident from these facts that in some re- 
gions of flow one must speak of the spectral content of the wake fluctua- 
tions in assessing the value of the Strouhal number and the forces exerted 
on the cylinder. The foregoing discussion of the Strouhal number raises 
two questions: firstly, is it possible to devise a universal Strouhal num- 
ber which would be applicable to all two and three-dimensional bluff 
bodies, and secondly, is it possible to theoretically predict the Strouhal 
number through, for example, the perturbation of the d'Alembert flow. An 
attempt along these lines has been made by Sacksteder Cl3 who found St = 
0.2028 for an infinitely large Reynolds number. The reasoning behind 
this prediction and its validity will not be discussed here further. 

The former question of the possibility of devising a universal Strouhal 
number has been of interest to various researchers. Table I, shown below, 
lists those definitions together with their originators. 



TABLE I 



St = 


Va 1 ue 


Body Type 


Origi nator 


f^h'/UK 


0.163 


f wo-d imens iona 1 bodies 


Roshko C 2 D 


f 1 /U 

V V 


0.181 


two-dimensional bodies 


Bearman C3l] 


f Jl /K‘U 
V w 


0.190 


ax i symmetric bodies 


Calvert 


f 9. /K‘U 
V w 


0.163 


two-dimensional bodies 


Simmons C53 



In Table I, h' represents the wake width as determined from Roshko's 
notched- hodograph theory Z2j; lateral spacing of the vortices as 
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determined by Kronauner's minimum drag criteria as reported by Bearman C3l] 

the wake width at a distance of the wake establishment region; and K = 

/l-C . where C . is the base pressure coefficient. All these definitions 
pb pb 

require either the measurement or the solution of one or more characteris- 
tics of the wake. Thus, their capability to predict a vortex shedding fre 
quency is quite limited. 

It is advantageous to note at an early stage those factors which in- 
fluence the Strouhal number. Firstly, small amplitude transverse oscilla- 
tion of cylinders or cables distinctly regularize the vortex shedding and 
hence the Strouhal number. For example, the experiments conducted at the 
National Physical Laboratory C63 have shown that the Strouhal number even 
in the drag crisis region remains perfectly constant and equal to about 
0.21. Secondly, the vortices do not peel off from a cylinder as perfectly 
two-dimensional line vortices. As it was first shown by Humphreys Ct 3, 
there Is a cell pattern separated by a correlation length. The longer the 
cylinder the larger the number of cells. The consequences of the forma- 
tion of such cells are that there are phase shifts in pressure along the 
cylinder and hence a reduction in the integrated transverse force. In 
other words, the transverse force calculated through the integration of 
the pressure distribution at a given section is not identical to the aver- 
age transverse force determined experimentally. All numerical analyses 
which model the complex three-dimensional flow as a two-dimensional flow 
must necessarily predict a transverse force or lift coefficient somewhat 
larger than those found by measurements. Other consequences of the cell 
formation concern the strength of the vortices. The vortices are not 
necessarily oriented in a plane normal to the flow and the measurements 
(flow visualization or electro-optical techniques) yield only that 
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component of vorticity which is normal to the plane of flow. Additional- 
ly, experimental results necessarily depend on the I ength-to-d i ameter 
ratio, end conditions, intensity and scale of turbulence in the ambient 
flow, the rigidity of the cylinder mounting, and of course on Reynolds 
number and the surface conditions of the cylinder. Most recent experi- 
ments (Shih and Hove CSD) have shown that roughness affects the base 
pressure in the same manner as the small transverse oscillations (beyond 
a threshold limit of about 0.10 diameters). Uniform roughness yields a 
perfectly uniform base pressure and removes the cell pattern. 

C. RESISTANCE IN STEADY FLOWS 
I . In-Line Forces 

The understanding of the fluid forces exerted on bluff bodies 
comes primarily from measurements. Of all bluff bodies, the circular cy- 
linder has been subjected to the largest number of investigations. Those 
experiments which constitute milestones in the determination of the magni- 
tude and characteristics of resistance are the works by Wiesel sberger C93, 

Fage and Johansen ClOj, and Roshko Cl lH. A plot of the drag coefficient, 

2 

Cp. defined by C_ = 2F/pLDU , is shown in Fig. I as a function of the 

U D '' 

Reynolds number defined by Re = UD/v. It is a well known fact that for 

4 5 

Reynolds numbers larger than about 10 and smaller than about 2x10 , the 
drag coefficient for a circular cylinder remains essentially constant at 

5 

a value of about 1.2. For Re > 2x10 , the shear layers emanating from 
the separation points become turbulent. Since a turbulent flow can ex- 
change larger momentum between the inner and outer flow (i.e., between the 
boundary I ayer and externa I flow), it can sustain a larger adverse pressure 
gradient and even reattach to the cylinder. The reattachment of the tur- 
bulent shear layers gives rise to the formation of so called separation 
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Drag coefficient versus Reynolds number for a smooth circular cylinder. 



bubbles and to the displacement of the attached flow further downstream 
from the first separation point. The reattached flow cannot continue to 
remain attached indefinitely and reseparates at an angle of about 130 de- 
grees from the front stagnation point. Thus, the transition of the flow 
in the shear layers from a laminar to a turbulent state results in a con- 
siderably smaller wake and hence in a considerably reduced drag coefficient. 
This particular flow regime is commonly called the drag crisis region or 
the critical region. As the Reynolds number increases, parts of the bound- 
ary layer adjacent to the first separation point become turbulent and the 
second separation point begins to move gradually upstream. This region is 
called the transcr i tica I region. At higher Reynolds numbers, the entire 
boundary layer becomes turbulent and the drag coefficient reaches a nearly 
constant value of about 0.6. There is considerable scatter in the data 
in the supercritical region because of the difficulties associated with 
flows at such high speeds. 

The present knowledge of the drag coefficient for a smooth cylin- 
der is limited to Reynolds numbers smaller than about 10^. There is great 
practical interest in the determination of force coefficients, i.e., drag 
coefficient, lift coefficient, and the Strouhal number for projects such 
as ocean thermal energy conversion (OTEC), large cooling towers, high rise 

structures, etc. A careful consideration of the various possible methods 

8 1 0 

to obtain data at Reynolds numbers in the order of 10 or 10 shows that 
there are practically insurmountable difficulties. The existing wind and 
water tunnels are not capable of yielding such high Reynolds numbers with- 
out the adverse effects of blockage and compressibility. The pop-up and 
drop-down experiments with large cylinders are constrained by cavitation, 
weight, and stability conditions. It might be possible in the future to 
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close the high Reynolds number gap by building large cryogenic wind tun- 
nels operating with freon or liquid helium. 

The drag coefficient depicted in Fig. I may be influenced by a 
number of factors such as roughness and turbulence. Numerous investiga- 
tors, notably Page and Warsap Cl2D, Achenbach Cl 3D, and Guven Cl 4], have 
clearly demonstrated that the larger the surface roughness, the sooner is 
the occurrence of the critical region and the higher the dip in the drag 
coefficient (see Fig. 2). Turbulence plays a similar role in the evolu- 
tion of the drag coefficients. It should be noted from Fig. 2 that the 
drag coefficient for a rough cylinder at supercritical Reynolds numbers is 
about twice that for a smooth cylinder. Furthermore, the experiments by 
Guven Cl4] and Szechenyi CiSD have shown that at sufficiently high 
Reynolds numbers there is a Reynolds number independent region. This ex- 
perimental fact has been termed the independence principle and used (Shih 

and Hove CSj) to extrapolate the drag coefficient data to Reynolds numbers 

0 

in the order of 1 0 . If the validity of the independence principle can 
be demonstrated unambiguously, it may obviate the need for difficult ex- 
periments in cryogenic tunnels. 

2. Transverse Forces 

Lift or transverse-force coefficient data are far less numerous and 

exhibit a great deal of scatter. Figure 3, which represents most of the 

existing lift coefficient data, shows that the lift coefficient, defined 

2 

by = (amplitude of lift force) /O. SpLDU , may vary from 0.2 to 1.0 in the 
subcritical Reynolds number range. In the critical range it is impossible 
to define a lift coefficient without resorting io its root-mean- square (rms) 
values. It is for this reason that some of the data appearing in Fig. 3 
are in terms of the rms values of the lift coefficients. For a more 
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Experimentally obtained values of lift coefficient versus Reynolds 
number for a circular cylinder Cl63. 



detailed discussion of the spectral content of lift in the critical region 
see Bub I i tz Cl 7D. 

The reasons associated with the large variations of the lift co- 
efficient are directly related to the mobility of the separation points, 
formation of vortex cells, lack of spanwise correlation, end effects, tur- 
bulence present in the ambient flow, etc. Since each experiment is carried 
out with a different I ength-to-d iameter ratio, blockage ratio, and end 
conditions, the number of vortex cells over the cylinder and the effect of 
the flow between the end of the cylinder and the tunnel wall give rise to 
widely different lift coefficients. 

The effect of roughness on the lift coefficient is to reduce its 
scatter and to increase its magnitude. As noted earlier, roughness dis- 
tinctly regularizes the vortex formation and the base pressure thereby 
eliminating the effect of phase shifts which cause the said scatter in 
lift. Data supporting the above experimental facts will not be reported 
here for sake of brevity. For additional details the reader is referred 
to the work by Szechenyi ClSD. 

D. RESISTANCE IN TIME-DEPENDENT FLOWS 

There are an infinite number of time dependent flows and it is impos- 
sible to discuss their characteristics in any unifying manner. One might 
obtain a clearer insight into such flows by considering only those cases 
which are relatively more manageable for numerical or laboratory experiments. 
Among such flows which have been subjected to extensive numerical analysis 
and experiments are the impulsively started flows, uniformly accelerating 
flows, and harmonically oscillating flows. Obviously, no distinction is 
made here between the said fluid flows about a body at rest and a body 
undergoing similar time dependent motions in a fluid otherwise at rest. 
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I . Impulsively Started Flows 



The impulsively started flow has been by far the one most investi- 
gated numerically through the use of the Navi er-Stokes equations. These 
investigations (Collins and Dennis CiSD, and Yang Cl9j) have been confined 
to the very early stages of the motion during which the separation points 
move from the rear stagnation point to about 109 degrees. An exception 
to the foregoing is the work reported by Thoman and Szewczyk C20D. They 
have solved the Nav i er-Stokes equations through the use of a suitable 
finite difference technique for various Reynolds numbers and found that 
the drag coefficient overshoots for a relative flow displacement (normalized 
time) of Ut/c = 5 and then gradually reaches steady state values. 

Only two separate investigations of the impulsively started flow 
have been reported. The first is that carried out by Schwabe C2lU upon 
PrandtI's suggestion and the other by Sarpkaya C22, 23j. Schwabe's experi- 
ments were conducted with a 9 cm diameter cylinder at a Reynolds number 
of about 600. His drag coefficient, calculated indirectly through the use 
of motion pictures, rapidly reaches a value of about 2.0 and remains con- 
stant thereafter. Schwabe's findings are not in conformity with either 
the numerical calculation or with the experimental results obtained by 
Sarpkaya. Figure 4 shows three sets of representative data at Reynolds 
numbers from 16000 to 35000 for the drag coefficient. Clearly, the very 
initial instants of the motion cannot be analyzed because of the problems 
associated with the generation of an impulsively started flow. At Ut/c 
slightly greater than zero, the force acting on the cylinder is almost 
entirely inertial. As the separation sets in at the rear stagnation point 
and rapidly moves upstream, the drag also increases rapidly and reaches 
a value of about 1.5 at Ut/c between 4.0 and 5.0. Then drops rapidly 
to the steady state value for smooth cylinders in steady uniform flows at 
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Ut/c larger than about 7.0. Detailed flow visualization experiments by 
Sarpkaya C22j have shown that the reason for the overshoot in Cj^ is the 
symmetric development of two vortices in the cylinder wake and the rapid 
accumulation of vorticity in each vortex. It must be noted that for such 
a symmetric vortex configuration, the only means by which the vorticity 
may be cancelled are the exchange of vorticity between the two vortices 
and the cancellation of vorticity with oppositely-signed vorticity, so 
called counter-vorticity, generated in the rear boundary layers. In other 
words, the rate of accumulation of vorticity in the vortices is consider- 
ably larger than that at any other time. Other vorticity cancellation 
mechanisms which come into play during the separation of the vortex from 
Its feeding sheet, as proposed by Gerrard C243, do not play any role dur- 
ing the growth of the symmetric vortices. Thus, the rapid accumulation 
of vorticity gives rise to the drag overshoot. As the asymmetry sets in, 
vortices shed a I ternati ng I y and the flow gradually becomes asymptotically 
steady. The possible mechanisms which cause the asymmetry in an otherwise 
symmetrical flow about a symmetrical body are discussed in greater detail 
later. Suffice it to note that one of the fundamental facts of fluid 
flows is that symmetric causes do not give rise to symmetric effects, as 
most eloquently philosophized by Birkhoff C25I1. 

Following the onset of asymmetry, the cylinder experiences a 
transverse force, as shown in Fig. 5, in terms of the lift coefficient and 
the normalized time Ut/c. The lift force fluctuates with the Strouhal fre- 
quency. The shedding of each vortex gives rise to a small fluctuation in 
Cp. The frequency of this fluctuation is twice the Strouhal frequency 
since two vortices (one from top and from bottom) are shed during a given 
cycle of lift osci I lations. 
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Fig. 4. Drag coefficient versus normalized time for impulsively-started flow about 
a circular cylinder C233. 
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Lift coefficient versus normalized time for impu I si veiy-started flow about 
circular cylinder C233. 



2. Other Time-Dependent Flows 



As noted earlier, the uniformly accelerating flows and harmonical- 
ly oscillating flows have been subjected to numerous analytical and experi- 
mental investigations. Among those notable are the investigations carried 
out by Iverson and Balent £26^, Kelm £21~], Hamilton C28j, Odar and Hamilton 
C293, Sarpkaya and Garri son C30H, and Sarpkaya C3l], The basic question 
has been and continues to be the formulation of a resistance equation which 
takes into consideration the effects of acceleration, velocity, and the 
history of the motion. All investigations have shown that resistance in 
time-dependent flows in general and in harmonically oscillating flows in 
particular cannot be considered as a juxtaposition of resistance in steady 
flows at corresponding instantaneous Reynolds numbers. Stokes, in his 
celebrated paper on the sinusoidal oscillation of a sphere in a fluid 
otherwise at rest, has shown that the fluid resistance is comprised of a 
linear velocity-dependent force and an acceleration-dependent inertial 

force. The coefficients for both forces depend on the Reynolds number and 

2 

the frequency parameter D /vT where T is the period of oscillation. Ob- 
viously, Stokes’ analysis is valid for only very small, as yet unspecified, 
Reynolds numbers. Inspired partly by Stokes' contribution and partly by 
practical considerations, Morison C32] proposed a resistance equation which 
may be written as 

F=lpApC^U|U| +pVC„^ (I) 

in which U represents the instantaneous velocity; A^, the projected area; 

V, the displaced volume; and C|^ and Cjy|, the drag and inertia coefficients 
respectively. For harmonically oscillating flow about a stationary cylin- 
der where the fluid velocity is U = cos 9 (here 0 = 2'irt/T) Eq. (I) 
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reduces to 



7 



- Cpicos 6|cos e + -JpiTp-Cj, sin 6 

M 



( 2 ) 



The essence of this equation is that the resistance in harmonic flow is 
assumed to be equal to a linear combination of a velocity-squared dependent 
drag force and an acceleration-dependent inertia force. In the midst of 
insufficiently clear insight into the kinematics of the complex flow, it 
has not yet been possible to demonstrate the reasons for the differences 
between the measured force and that calculated from Morison’s equation for 
U|^T/D values in the vicinity of 15. Additional experiments and numerical 
analysis, possibly through the use of the discrete vortex model, might shed 
further light on the subject and into the overall understanding of resis- 
tance in time-dependent flows. It is evident from the discussion of both 
the sj'eady and unsteady flows that analytical or numerical methods are 
urgently needed with which relatively inexpensive investigations can be 
conducted to clarify the complex wake-boundary- layer interaction mechanism. 
Of all the possible flows which might be considered as a candidate for such 
an analysis, the impulsively started flow is certainly the one most manage- 
able and most fundamental to the understanding of a I I other flows. 



E. METHODS OF NUMERICAL ANALYSIS 

With the advent of high speed computers, numerical methods began to 
play a greater role in fluid dynamics. Currently four major methods of 
analysis are in use: finite difference, finite element, marker-cell tech- 

nique, and the discrete vortex method. 

Finite difference methods have been used for the solution of flow 
about a circular cylinder for relatively small Reynolds numbers (Re = 100) 
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either through the use of the Nav ier-Stokes equations as they are commonly 
written or through the more frequently used vorticity equations. Ingham 
C33] calculated the steady flow past a circular cylinder at Re = 100. 

Thoman and Szewczyk C20j, as noted earlier, carried out a similar numeri- 
cal analysis of flow over stationary and rotating cylinders at Reynolds 
numbers as large as 10^. There are a number of such examples in the 
literature but they are not cited here because of their lack of relevance 
to the subject under consideration. 

The finite element method, which came into prominence first in the 
area of structural analysis and then permeated through heat transfer and 
fluid mechanics, has been primari ly used for unseparated flows about 
bluff bodies. As far as fluid mechanics is concerned, the outstanding ex- 
amples of such applications may be found in two recent volumes of work 
dealing with finite elements in fluid mechanics C34H. Bretanow and Ecer 
C35D made an attempt to calculate the separated flow about a circular 
cylinder at small Reynolds numbers but the results have been inconclusive. 
Applications of the finite element method to flows about foils where the 
flow remains attached have yielded results in conformity with experiments 
(Bretanow and Ecer C35l]). 

The marker-cell technique developed by the researchers at the Los 
Alamos Scientific Laboratory during the past 20 years uses the finite 
difference form of the Nav ier-Stokes equations and then adjusts the pres- 
sure in each cell to insure that the equation of continuity is satisfied. 
Numerous ingenious techniques have been developed by Fromm and Harlow 
C36l] and Nichols and Hirt C37l] for the efficient adjustment of pressure 
in each cell. Such techniques will not be dealt with here. The researchers 
of Los Alamos Scientific Laboratory have applied the marker-cel I technique 
to the flow past cylinders, spheres, rectangular blocks and to various other 
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free-surface flows. These applications dealt with very low Reynolds num- 
bers (in the order of 100). Even though the technique has considerably 
increased the understanding of the kinematics of the wake, the precise 
nature of the time-dependent force acting on any one of the bodies cited 
above remains obscure. 

The discrete vortex model, which is included as one of the four major 
models, is the subject of extensive discussion in the next chapter. There- 
fore, the remainder of this section is confined to a discussion of the suc- 
cesses and shortcomings of the first three methods. 

All numerical models are by their very nature approximate. They satis- 
fy certain critical conditions only at discrete time steps. Secondly, they 
require mesh or element sizes which vary from one region of flow to another 
to insure that the characteristic gradients of the flow (such as velocity 
or pressure) do not significantly or abruptly change over any one element. 
For example, the finite difference method when applied to the viscous flow 
over a circular cylinder requires that the size of the elements on and near 
the cylinder be considerably smaller than the boundary layer thickness. 
Otherwise, the evolution of the boundary layer, the positions of the sepa- 
ration points, and the evolution of the wake cannot be accurately calcu- 
lated. Since the boundary layer thickness for laminar flow is proportional 
- 1/2 

to Re , one realizes that only for relatively small Reynolds numbers can 
a sufficiently large boundary layer thickness and a correspondingly large 
mesh size be obtained. This fundamental difficulty has prevented the 
finite difference, finite element, and the marker-cel I technique from making 
significant progress into the calculation of steady or time dependent flows 
about bluff bodies at Reynolds numbers of practical interest. None of the 
three techniques cited above has been applied to the solution of flow prob- 
lems associated with oscillating bodies. The reasons for this are rather 
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self-evident since the solution requires not only a moving grid but also 
a grid whose mesh size must be varied with time to minimize the effect of 
sharp velocity and pressure gradients. Finally, it should be noted that 
all three numerical techniques present "stiffer" solutions and introduce 
an artificial numerical viscosity whose influence cannot be readily 
assessed. Not to be entirely critical of the said methods, for there are 
correspond i ng shortcomings of the discrete vortex model, one must note 
that in the wide spectrum of Reynolds numbers there will always be some 
interest in low Reynolds number flows. Thus, such numerical techniques 
will continue to play some role in separated flows where the discrete 
vortex model is not advantageously suited. It therefore appears that no 
one model will stand out as superior and universally applicable and that 
all methods will continue to complement each other. 
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11. FUNDAMENTALS OF THE DISCRETE VORTEX MODEL 



A. INTRODUCTION 

The discrete vortex model (DVM) is a potential flow representation of 
the shear layers emanating from the separation points. Any separated vis- 
cous flow is completely specified if the spatial and temporal distribution 
of vorticity in the entire flow field is known. For reasons described in 
the previous chapter, it has not yet been possible to achieve this objec- 
tive through any exact or approximate methods. It has been observed and 
experimentally verified by Fage and Johansen ClOj that the shear layers 
at sufficiently high Reynolds numbers are quite thin and the vorticity is 
confined into regions of tightly spiralled vortex sheets. This fact alone 
invites one to divide the shear layer into a number of small segments and 
concentrate the vorticity of each segment into a line vortex situated at 
the geometric center of the segment. This idealization permits one to 
use the powerful complex function theory to determine the kinematics and 
dynamics of the flow. Evidently, the validity of the results can only 
be judged in light of experimental facts. 

In the following, only the gross features of the method are described 
and the details are presented following an historical resume. 

At a given time the separated shear layers are represented by vortices 
situated at the center of each segment. The potential flow theory permits 
one to calculate the velocity of each vortex induced at its center by all 
other vortices and the ambient flow. Then the vortices are convected for 
a short time interval At through the use of a suitable convection scheme 
to new positions. It is then necessary to introduce new vorticity into 
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the flow at or near the separation points to account for the vorticity 
generated in the boundary layer. The introduction of new vortices (nascent 
vortices) represents an addition to each set of vortex arrays. Then the 
process is repeated. 

In spite of its apparent simplicity, there are a number of problems 
which must be resolved in carrying out the said calculations: (i) the vor- 
ticity introduced at each time step must be related to the characteristics 
of the boundary layers upstream of the separation points; (ii) one must 
insure that the point vortices continue to represent the vorticity at each 
segment in accordance with the fundamentals of hydrodynamics; (iii) allow- 
ance must be made for the interaction of oppositely-signed vorticity and 
the decay of vortices; (iv) methods must be devised to limit the number of 
the point vortices to minimize the computation time without affecting the 
kinematics and dynamics of the flow; and (v) procedures must be established 
for the flow to retain its inherent flexibility without introducing 
artificial stiffening effects (e.g., to allow for the mobility of stagna- 
tion and separation points). 

The foregoing conveys the idea that the discrete vortex model is and 
will continue to rely on a strong interaction between potential and vis- 
cous flow theories. Neither is capable of providing a solution alone but 
together they can make inroads towards describing the basic features of 
the separated flow provided that sufficient ingenuity is incorporated into 
the said interaction. Lastly, it must be noted that all the numerical 
methods described above depend on the availability of a high-speed digital 
computer. 

B. HISTORICAL DEVELOPMENT 

The DVM has been in existence for the past 46 years since its incep- 
tion by Rosenhead C38] in 1932. Rosenhead investigated the so-called 
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"Helmholtz instability" of a surface of discontinuity, separating two 
fluids of equal speed and density but moving in opposite directions. His 
primary objective was to give the shear layer between the two fluids an 
initial disturbance by situating the shear layer along the curve y = 0. I 
sin (27 tx/X) and then observe the evolution of the sheet with time. To ac- 
complish this, Rosenhead represented the shear layer by point vortices 
placed along the sheet at intervals such that their projections along the 
x-axis were equ i-d i stant. Each vortex was assigned a strength of UX/6n 
where n is the number of vortices per wave I ength, X, and U the speed of 
the two fluids. Keeping the vortex strengths constant, Rosenhead followed 
the deformation of the sheet by convecting the vortices with a simple 
Eulerian integration of t ^;|^At where T'* is the position vector 

of the j-th vortex; q'^, the computed velocity vector; and At, the time 
interval. Using n = 2, 4, 8, and 12, Rosenhead observed that in each case 
the sheet rolled up smoothly into concentrated clusters of vortices posi- 
tioned at intervals of X along the x-axis. Rosenhead’s results, which were 
computed using a desktop calculator, are shown in Fig. 6. The tendency of 
a distribution of point vortices to roll up in a manner consistent with 
observations is clearly evident, 

Rosenhead's results remained undisputed until 1959 when the computa- 
tional power of the digital computer was applied to the analysis. Birkhoff 
and Fisher C39l] repeated Rosenhead’s analysis using a Runge-Kutta integra- 
tion scheme, 22 vortices per wavelength, and smaller time intervals. In- 
stead of improved results they observed that the vortices took irregular 
paths and no longer rolled up smoothly. Birkhoff and Fisher concluded that, 
given sufficient time, an originally organized array of equ i-strength point 
vortices would eventually become random and that the validity of represent- 
ing a shear layer by point vortices is questionable. 
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Fig. 6. Example of Rosenhead's 
calculations C38]. 
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Shortly thereafter, Hama and Burke undertook a similar computer 

study and verified the observations of Birkhoff and Fisher. However, Hama 
and Burke noted a conceptual error in Rosenhead's original method of plac- 
ing the discrete vortices along the shear layer. Specifically, if the 
initial vortex sheet was to have a constant vorticity per unit distance 
along its length and be represented by equ i-strength vortices, then the 
vortices should be placed at equal distances measured along the sheet vice 
at equal x increments. Using this arrangement Hama and Burke repeated 
their calculations and found that the onset of irregular vortex motion was 
significantly delayed. Eventually, however, a random pattern did result. 
Thus, at this point it was hypothesized that the onset of instability 
could be delayed although the questions of why it occurred and how it 
could be controlled were still a mystery. 

Without fully realizing the consequences of sheet instability (the 
result of which is discussed later), the DVM continued in its evolution 
relying on the assumption that meaningful results could be obtained up to 
the onset of Irregular vortex motions. In 1962 Abernathy and Kronauer C4lD 
applied the method to investigate the motion of two shear layers of oppo- 
sitely-signed vorticity separated by a distance h^ and each given an 
initial perturbation as in Rosenhead's work. They observed that the vor- 
tices accumulated into clusters resembling the pattern of a vortex street. 
An example of their results is shown in Fig. 7. 

Encouraged by the demonstrated ability of the DVM to represent a vor- 
tex street, the next logical extension of the method was to analyze the 
flow past a bluff body. However, the inclusion of a solid body in the 
flow field necessitated new considerations. Specifically, the previous 
analyses had begun with a discretized sheet having a preassigned initial 
position and vorticity distribution along its length (usually taken as 
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constant). The application of the DVM to flow past a bluff body required 
that the discretized shear layers be generated continuously by adding 
additional point vortices into the flow at regular time intervals and at 
locations representative of where the shear layers actually emanate from 
the body. The requirement for some method to account for vorticity genera- 
tion and introduction in an otherwise potential flow model posed a signi- 
ficant problem. Compounding the difficulty was the time dependent nature 
of both the generation rate and its point of introduction. 

The first DVM analysis to include a bluff body was performed by 
Gerrard C42H for a circular cylinder. Although a circular cylinder was of 
great practical interest it was also the most difficult one to analyze 
because of the unsteady vorticity generation and mobile separation points. 
Gerrard was forced to adopt a rather arbitrary method of introducing new 
vortices which severely influenced the results of his analysis, although 
it still represented a noteworthy application of the DVM. Noting the dif- 
ficulties associated with the mobile separation points and the vorticity 
introduction, Sarpkaya C433 introduced a mobile separation region repre- 
sented by an infinite number of vortex sheets of vanishingly small vorticity 
connected to the nascent vortex. This analysis was restricted to impul- 
sively started symmetric flow about a circular cylinder. In the ensuing 
years many i nvest igators bypassed the difficulty introduced by the mobility 
of the separation points and dealt with bodies with sharp edges (e.g., 
inclined plates, blunt-based bodies, wedges, etc.) where the separation 
points are fixed. Despite the growing number of applications of the DVM, 
a number of disturbing difficulties remained unresolved and a full inter- 
action between the boundary layer, separation points, and the wake had 
not been established. These difficulties required the use of a number of 
nondisposable parameters. Consequently, the proper balance between 
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heuristic reasoning and computational formulation had to evolve through 
the recognition of the reasons underlying the basic difficulties. 

Briefly, the most important of these difficulties stemmed from the 
proximity of the vortices to each other and to the solid boundary, and 
from the necessity of representing the continuous process of vorticity 
generation at the singular separation points by a discrete process. The 
occasional proximity of two point vortices resulted in large mutually- 
induced velocities and hence the shooting of the two vortices rapidly 
away from the flow field. This was particularly true when the vortices 
ceased to be located at the geometric center of the sheet segment which 
they were supposed to represent. Furthermore, the spiralling of the sheets 
brought vortices on two adjacent spirals close to each other and gave rise 
to a traveling and growing instability. Thus, it was not possible to 
correctly represent the core of tightly spiralled sheets. The proximity 
of the vortices to the boundary and hence to their images inside the 
body resulted not only in similar instabilities but also in unrealistic 
local pressures. This, in turn, resulted in unnatural fluctuations in 
drag and lift forces. 

Faced with the problems just cited and the desire to enhance the 
applicability of the method, various investigators devised several numeri- 
cal cures. Chorin and Bernard C44] introduced a cut-off distance from 
the center of each vortex beyond which the vortices were to behave like 
potential vortices. For radial distances smaller than the cut-off distance 
the vortices were assumed to exhibit a solid body rotation. Their calcu- 
lations did not, however, result In an improvement of the analysis of 
separated flow about a circular cylinder. Moore C45H, in dealing with 
two vortex sheets emanating from a plane lifting surface, introduced the 
method of amalgamation. When a vortex sheet was spiralled tightly enough. 
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the last vortex on the sheet was connected to a core vortex. At each time 



step the last vortex was amalgamated into the core vortex increasing its 
strength. This procedure partly avoided the difficulties stemming from 
the proximity problems of the adjacent coils of the spiralling sheet. 

Both investigations dealt with the devising of methods to overcome the 
instability problems rather than with the understanding of the real cause 
of the d i f f i cu I ty. 

The basic question of why an initially equidistant set of vortices 
should after several convections come in close proximity of each other was 
resolved by Fink and Soh C46H who demonstrated that the primary source of 
the difficulty was in failing to rediscretize the sheet at each time in- 
terval of the calculation. In other words, it is necessary to insure that 
the point vortices are always located at the geometric center of the seg- 
ment of the sheet which they are supposed to represent. The complex con- 
jugate velocity q(z.) at a point z. due to a single vortex of strength F. 

J J ^ 

located at is given by* 



q(z . ) 
J 



2i7T z . - z, 



(3) 



For a series of N discrete vortices on a vortex sheet, the same velocity 
at a point z. becomes C463 
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(4) 



Mathematical details of the relevant potential flow theory will be 
presented in Chp. III. Here no attempt is made to introduce the precise 
definitions of the variables appearing. 
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Equation (4) resembles Eq. (3) for the induced velocity field due to point 

vortices except for the logarithmic term. The point z. lies within the 

J 

segment (Sj^j^ 2 » without necessarily bisecting it at all times, and 

s measures the distance along the sheet. The consequences of Eq. (4) are 
that: (i) if the equivalent point vortex is not placed at the midpoint of 
its segment through rediscretization of the sheet at each time interval, 
then the logarithmic term does not vanish and the computational error in- 
creases depending on the number of vortices, the time interval used, and 
the total time of computation; (ii) the vortices, which initially bisect 
the segment which they are to represent, do not continue to do so at the 
succeeding time intervals; (iii) the use of finite vortex cores, amalgama- 
tion of vortices at the center of the spiral, or other techniques only 
delay or minimize the accumulation of the errors resulting from the loga- 
rithmic term in an amount related to the distance between z. and the center 

J 

of the segment; and that (iv) the growth of the computational error may be 
significantly reduced by placing each discrete vortex at the midpoint of 
its segment, i.e., by placing the vortex at Zj = 0.5(Zj_|^|y2+Zj_ S'f’ each 
time interval. Only through such a procedure can one make the logarithmic 
term vanish. 

The calculations are then carried out at each time step by represent- 
ing the vorticity density by an entirely new set of equidistant vortices 
whose strengths are adjusted to give a good representation of that density. 

The foregoing procedure does not resolve all of the computational 
errors particularly in regions where the radius of curvature of the sheet 
is small. Fink and Soh applied this technique to the early stages of 
relatively simple flows wifh fixed separation points. Their examples did 
not deal with the alternate shedding of vortices. 
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Even though the method of rediscretization essentially resolved the 
problems stemming from the proximity effects, two fundamental links between 
the potential flow model and the behavior of real fluids remained unre- 
solved. The first of these is the relationship between the vorticity to 
be introduced at each time step and the boundary layer which generates it. 
The second is the decay of vortices by various mechanisms. The first 
question was initially considered by Deffenbaugh and Marshall C473 who in- 
vestigated the Impulsively started flow about a circular cylinder (without 
recourse to red i scret i zat ion ) . They used Pohihausen's boundary- I ayer ap- 
proximation method to determine the position of the separation point and 
the vorticity to be introduced at each time interval. Their analysis pre- 
dicted a separation angle of 6^ = 64 degrees, which is considerably smaller 
than the commonly accepted experimental value of 9^ = 81 degrees. - 

The question concerning the decay of vortices has remained practically 
obscure. Although several researchers employed methods to cancel opposite- 
ly-signed vortices whenever they came within a prespecified distance of 
each other, no systematic investigation of the effect of decay of vorticity 
was conducted. Often recourse was made to PrandtI's suggestion that indi- 
vidual vortices in a vortex street retain only about 60 percent of the 
vorticity generated in the boundary layer during a period required to shed 
a single vortex. Some investigators, in efforts to obtain results in 
partial conformity with PrandtI's suggestion, automatically applied a 40 
percent reduction to the strength of each nascent vortex. Seldom was men- 
tion made of the experimental ly observed fact that the rate of decay of 
vortex strengths varies with distance behind the cylinder. 

The foregoing represents the state-of-the-art prior to the beginning 
of the present investigation. A relatively more detailed account of the 
historical evolution of the DVM is found in Clements and Maul I C48]] and 
Fink and Soh [461]. 
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III. THEORETICAL AND PHENOMONOLOGICAL FOUNDATIONS 
OF THE DISCRETE VORTEX MODEL 



A. INTRODUCTION 

This chapter deals with the characteristics of inviscid and viscous 
flows about a circular cylinder from the forward stagnation point to the 
vortices in the far wake. Mathematical principles, boundary layer con- 
cepts, and relevant hypotheses are introduced in preparation for the dis- 
cussion of the details of the discrete vortex model developed in the 
present study. 

B. FLOW KINEMATICS 

The complex potential for a two-dimensional inviscid flow is given by 

w(z) = <)) + iip 



where (p and ip represent, respectively, the velocity potential and the 
stream function. The velocity components which satisfy the Cauchy-Ri emman 
conditions are given by 



9x 9y 



and 



= fli 

9y 9x 



In terms of the complex potential, one has 



(5) 



dw(z) 

dz 



u + iv 



( 6 ) 



The flow about a circular cylinder is generated by combining a uniform 
flow with a doublet at the center of the cylinder (see e.g.. Mi I ne-Thomson 
C49j). This yields 

w(z) = w = - U(z + c^/z) (7) 
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in which -Uz represents the uniform flow and the remainder the doublet. 
The introduction of a vortex into the flow field exterior to the cylinder 
requires, by virtue of the circle theorem, the introduction of an image 
vortex interior to the circle. Thus, the flow about a circular cylinder 
with N number of real vortices is given by 



w 



-U(z + — ) + 



N 

sE 

n= I 



r {Jln( 2 -z ) - Zn(.z - ■=— )} 
n n 7 



( 8 ) 



It should be noted that a real vortex located at z = z^ has an image at 

z = c /z^ and at z = 0. The image at the center of the circle, which 

gives rise to terms like FJlnCz), must be excluded due to the condition 

that the vortices shed from the cylinder leave a circulation opposite to 

their own on the cylinder. Clearly, the absence or presence of an image 

vortex at the center of the circle does not alter the boundary condition 

0({)/8n| = 0) since the streamlines of a vortex are circles. No re- 

|z|=c 

striction has been imposed regarding the time dependence of the velocity 
U other than that it be uniform. 

Substituting Eq. (8) into Eq. (6), one obtains the velocity at any 
point z as 



clw ,,,, c i 

-U+,V =_ = -U(|-^)+2y 



V r {— 

Z-j n z-2 



n=l 



c 

z - ■=- 
z_ 



(9) 



Equation (9) may also be used to evaluate the convective velocity of a vor- 
tex at z = z provided that the singular term F / iz-z) is omitted by virtue 
n ^ ^ n n ' 

of the condition that a vortex cannot impart a velocity on itself. 
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It is advantageous to work with dimensionless parameters by introduc- 
ing the following change of variables: 

z' = z/c, u' + iv' = (u+iv)/u, r' = r/uc (lo) 

w'(z) = w(z)/Uc t’ = Ut/c 

Substituting into Eqs. (8) and (9), letting U = 1.0, c = 1.0, and dropping 
the primes for sake of simplicity, one has 



N 

wCz) = -(z+l) + 4-Y^ r {5-n(z-z ) - £n(z-l/z )} 
z 2n n n n 

n= I 
and 



-u+i V 



2 




r_{i 



n z-z 



— L-) 

Z-I/-Zn 



(II) 



( 12 ) 



where u, v, F and z are now the normalized variables and w the normalized 
complex potential function. Hereafter, all measures of distance, time, 
velocity and circulation are given in terms of the above normalized 
variables. Any exceptions are clearly specified. 

Equation (12) is used extensively In convecting the vortices and in 
calculating the lift and drag forces. In fact, about 90 percent of the 
computation time of the present model is consumed in the evaluation of 
complex velocities. This points out the necessity of keeping the number 
of discrete vortices as small as possible while maintaining a sufficiently 
small time interval to insure that the character i st ics of fiow are faith- 
fully reproduced. The selection of the time interval, like the selection 
of the mesh size, is one of the subjective items of the method and 
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requires a careful parametric analysis. This is addressed to i n a later 
section. 



In general, the calculation of velocities and the convection of vor- 
tices are straightforward and noniterative. Under certain circumstances 
vortex-vortex proximity between adjacent spirals and the vortex-boundary 
proximity may require special attention. The resolution of such secondary 
issues is deferred to another section since they are not of any major con- 
sequence either theoretically or practically. 

C. FORMULATION OF THE RESISTANCE EQUATIONS 

The forces acting on a stationary circular cylinder of unit length by 
a time-dependent flow containing N vortices of strengths and locations 
can be calculated from the normalized form of the generalized Blasius 
theorem (Mi I ne-Thomson C49]) to give 



Inserting Eq. (II) into Eq. (13) and evaluating the residues, Sarpkaya 
CSOD has shown that 




(13) 



where and represent the drag and lift coefficients respectively. 



N 




(14) 



in which u^ and v^ represent the velocity components of the n-th real vor- 



tex and u^' and * , the velocity components of the image of the n-th vortex. 



Equation (14) may be written in a more compact and computationally 



more suitable form as 



N 




(15) 



n=l 
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where the dot denotes derivative with respect to time such that 



2 =q=u+iv (16a) 

n n n 



and 



• I I < . • < 

2 = q = u + IV 

n ^n n n 



(16b) 



The complex velocities q^ and q^ may easily be shown to be related by 



q 



i 

n 




(17a) 



or by 



+ iv 



-u 

n 



+ iv 

n 




(17b) 



The norma I i 2ed force given by Eqs. (14) or (15) can also be evaluated in- 
directly through the use of surface pressures and the unsteady form of 
Bernoulli's equation. In either case the results are identical and the 
use of one method or the other is simply one of personal choice and pro- 
gramming efficiency. 



D. GENERATION AND INTRODUCTION OF VORTICITY 

The vorticity per unit length of a thin shear layer separating two 
streams with velocities Uj and U 2 is equal to the difference of the velo- 
cities U|-U 2 . This vorticity is transported with the average velocity of 
the two streams, i.e., (Uj + U2)/2. This is a continuous process and the 
amount of vorticity transported during a time At becomes 

AT = (U|-U2)(U|+U2) At/2 = (Uj^-U2^)/2- At (|8) 
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The discrete vortex model represents the continuous generation and flux 
of vorticity into the outer flow by introducing point vortices at discrete 
time intervals. The strength and position of introduction of each vortex 
must depend on the interaction of the unsteady vorticity distribution in 
the wake and on the appropriate boundary- I ayer considerations necessary to 
determine where the separation occurs. Those applications of the DVM 
which consider the separation points fixed ignore not only the said inter- 
action but also artificially constrain the evolution of the flow. 

In the present analysis the instantaneous position of the separation 
points is determined through the use of an approximate boundary (ayer 
analysis. At a given instant the velocity distribution on the forebody 
of the cylinder is calculated and the corresponding pressure distribution 
is assumed to be impressed on the boundary layer. This fundamental feature 
of the boundary layer theory together with the integral momentum equation 
of von Karman permits the use of a number of separation prediction methods. 
Among the most notable of such methods are those given by Poh I hausen CSiX 
Stratford C52], Thwaites C53H, Curie and Skan C54], Timman C55D, and 
Loitsidnski E563. am of these methods exhibit varying degrees of sensi- 
tivity to the external velocity distribution and may not all yield identical 
or uniformly consistent results. 

It is important to note that the combination of potential and viscous 
flow methods cannot establish a one-to-one relationship between the Reynolds 
number and the flow character istics. It is only the kind of separation 
criteria used (laminar or turbulent) that determines whether the flow is 
subcritical, critical, or turbulent. The fact that the use of the DVM be- 
comes increasingly more accurate with decreasing shear layer thickness 
makes one realize that the results obtained with a laminar separation 
criteria correspond to moderate to high subcritical Reynolds numbers. 
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The i n teg ra I -momentum methods such as that of Poh I hausen require that 
the velocity outside the boundary layer be given and that the separation 
points be either fixed or fluctuate with small amplitudes and velocities 
about a mean position. In other words, the changes in the outer flow and 
the separation points should be such that the flow may be treated as 
quasi-steady. 

The outer velocity distribution may be calculated either numerically 
through the use of an interactive scheme between the wake and the boundary 
layer or through the use of experimentally determined pressure distribution. 
Hiemenz C57H used the measured pressure distribution about a circular 
cylinder and approximated the velocity distribution at the edge of the 
boundary layer by 

U(?) = 1.8348? - 0.27422?^ - 0.04782?^ (19) 

where ? represents the arc length measured along the cylinder from the 
stagnation point. Hiemenz's solution indicates that separation occurs at 
an angle 0^ = 82 degrees, measured from the stagnation point. This cele- 
brated result has often been taken as proof that in subcritical flows 
separation occurs at 6^ = 82 degrees. The fact has not been realized that 
the entire velocity distribution and in particular that portion between 
the maximum velocity and the separation point (i.e., from 6 - 70 to 0 - 
80 degrees) dramatically affect the position of the separation point. In 
fact, Meksyn C583 has shown that the higher order derivatives of Hiemenz's 
velocity profile do not satisfy certain fundamental criteria regarding 
the separation point. At present there is no completely satisfactory velo- 
city distribution for the forebody of a circular cylinder either in sub- 
critical or in the supercritical regimes. Additional discussion of this 
subject is found in Meksyn C58j. 
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Numerous experiments have yielded laminar separation points ranging 
from 0^ = 70 degrees to 6^ = 90 degrees. The commonly accepted value is 
about 80 degrees. Page and Johansen's ClOU experiments, and others since 
then, have shown that the separation point is in reality a separation zone 
which extends over a finite region and that spanwise and chordwise fluctua- 
tion of this region do not permit a precise determination of the mean 
position, amplitude, and period of oscillation of the separation point. 

The purpose of the foregoing brief discussion is to point out both the 
theoretical and experimental difficulties encountered in the determination 
of the separation point and the inadequacies of the appropriate laminar 
boundary layer methods. The prediction of the separation of turbulent 
boundary layers is even more uncertain (Stratford C52j, Takada C59]]) and 
are not discussed herein. Instead, attention is focused on Pohihausen's 
method for time-dependent boundary layers. The former is more appropriate 
to the flow conditions where the motion of the separation point and the 
velocity at the edge of the boundary layer may be treated as quas i -steady. 
Schuh's C603 method is more appropriate for the initial instants of motion 
where the separation points move rapidly from 6^ “ degrees to about 
0g = NO degrees. The motion of the separation point commences after a 
relative fluid displacement of (Schlichting CSlD) 

+ 'in-r473?T- “-^51 (20) 

The reason for this is that in impulsively started flows about rounded 
bluff bodies the separation does not begin until the fluid moves some dis- 
tance past the body. The only exception to that is the flow about bodies 
with sharp corners in which case the separation begins at t = 0.0. 

For the case of impulsively started flow about a cylinder, the time 
in the boundary layer calculations is measured from the instant of the 
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inception of separation. Thus, a relative time of t = 0.351 must be added 
to that used in the boundary layer calculations in order to match the 
starting times. 

Pohihausen’s method Is well described in the literature (Schllchting 
C6ll]) and only the highlights of the method are presented herein. Assuming 
a fourth order polynomial for the velocity distribution within the boundary 
layer, Pohihausen has shown that von Karman’s integral momentum equation 
may be reduced to 



dZ _ F(K) 

dC U(C) 



and 



K = 



7 dU(C) 
^ dC 



( 21 ) 



in which 



K = ( 



37 I 
315 ■ 945 



A - 



9072 



2 2 
A 



(22) 



with 

. _ 6^ dU(C) 

V dC 



(23) 



Other variables appearing in Eqs. (21), (22) and (23) are 



U(C): +he normalized velocity on the cylinder 

the arc length measured along the cylinder from the 
stagnation point 

6: the boundary layer thickness normalized by the cylinder 

radius 

V: viscosity normalized by the ambient velocity and the 

cylinder radius (sometimes referred to as the unit 
Reynolds number). 

The values of K and F(K) are tabulated in terms of. A (Schllchting CSlD). 
Similar tables are provided for other boundary layer approximations cited 
previous I y. 

Integration of the coupled Eqs. (21) begins at the front stagnation 
point and proceeds along both sides of the cylinder independently until 
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the separation criteria (A = -12) is satisfied for each layer. This re- 
sults in two distinct separation points which in general are not at 
identical angular distances from the mobile stagnation point. 

As noted earlier, the foregoing method is applicable only to the rela- 
tively established flow where the velocity at the edge of the boundary 
layer and the separation point do not significantly vary with time. During 
the early stages of motion, however, the time dependence of the said para- 
meters is very strong and must be treated through the use of the time 
dependent integral boundary-layer equations. Schuh C603 has extended 
Poh I hausen's method for the general calculation of the unsteady boundary 
layer. The details of his method are not presented here. The results 
obtained by Schuh for an impulsively started flow about a circular cylin- 
der are shown in Fig. 8. These results, based on the classical potential 
flow velocity distribution about a circle, show that the separation point 
rapidly approaches a value of about 6^ = 109 degrees. Schuh 's method is 
not capable of incorporating the growth of the symmetric wake and the 
subsequent motion of the separation points and thus is not valid at large 
times. Thoman and Szewczyk H20H employed the finite difference form of 
the full Nav ier-Stokes equations and obtained results identical to those 
shown in Fig. 8 for the initial instants of flow. They have further shown 
that the separation po i nt once havi ng reached an angular position of about 
110 degrees moves gradually to the final value of the separation angle 
(about 80 degrees). Partly to confirm this conclusion and partly to assess 
the significance of time dependence on the instantaneous position of the 
separation points, two sets of calculations are performed. Both use the 
time-dependent velocity distribution available through the model and both 
begin with an initial separation point of 6^ = 109 degrees. The first cal- 
culation uses Pohl hausen's steady-state separation criteria with the in- 



53 




Fig. 8. Separation angle versus time [603. 
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stantaneous velocity distribution. The second calculation uses the unsteady 
i n teg ra I -momentum equations given by Schuh C60H, including those terms 
associated with the unsteady velocity distribution. The results obtained 
from both sets of calculations are shown in Fig. 9. Evidently the two cal- 
culations yield almost identical results, pointing out the fact that the 
separation point moves rather slowly and that the contribution of the time 
dependent terms in the unsteady form of the equations are negligible. It 
is because of this reason that all calculations are based on Pohihausen’s 
separation criteria assuming the flow to be quas i -steady. This assumption 
improves as one approaches steady state. Figure 9 indicates that the 
angular speed of the separation point is about 10 degrees per second. The 
same speed about the mean position of the separation point is about 2 de- 
grees per second. (As is shown later, the separation point oscillates 
about a mean position of 0^ = 77 degrees with an amplitude of 3 degrees 
and a period of about 10.) Based on the foregoing, an amplitude of oscil- 
lation as much as 15 degrees can be performed without the inclusion of 
time-dependent terms. 

In the present model the calculations begin at the time of t = 1.0 
when the separation point is at 6^ = 109 degrees. It should be noted 
that a stra Ightforward application of Pohihausen's method with a potential 
flow velocity distribution (U = 2 sin 6) yields a separation angle of 
0^ = 109.5 degrees. 

Once the position of the separation points is determined, the rate of 
generation of vorticity is easily determined by writing 



6 




(24) 



o 
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Fig. 9 . 



Comparison of steady and unsteady 
separation point calculations. 
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where 6 is the boundary layer thickness; u, the velocity within the bound- 
ary layer; and y, the distance along the outward normal. The contribution 
of 8v/8x is excluded from the vorticity expression in conformity with the 
approximations of the boundary layer theory. 

The integration of Eq. (24) to obtain the circulation flux at a separa- 
tion poi nt yields 




dt 2 



(25) 



where is the tangential velocity at the predicted separation point. 
Ordinarily the rate of change of circulation should be written as 




dt 2 



(26) 



where Uj and U 2 represent respectively the velocities at the outer and 
inner edges of the separating shear layer. Extensive experiments with 
various bluff bodies by Page and Johansen ClOU have shown that is about 

dr 2 

5 percent of Uj and that is correctly represented by /2. This is a 

universally accepted experimental and theoretical fact. 

dr 

In numerical calculations ^ is replaced by AF/At and the strength of 
the nascent vortex is calculated from 
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nv 
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U At 
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(27) 



The placement of nascent vortices into the flow to adequately represent 
the separating shear layers has been a source of difficulty. This is pri- 
marily because the discretization of a continuous process near a singular 
region (separation point) is an exceedingly complex problem. In the present 
analysis, no new procedures are needed for the determination of the strength 
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and position of the nascent vortex because the separation point is obtained 

from the boundary layer analysis. By virtue of its definition, the 

separation point must be the point at which both the velocity and the 

2 

shear stress on the cy I i nder are zero . Thus, the nascent vortex of known 
strength must be placed on a radial line passing through the separation 
point at a distance m, i.e., at = (l+e) exp Ci(‘fT-6^)D, to render the 

tangential velocity on the cylinder at the separation point equal to zero. 
This yields a unique value for e as 






(28) 



The use of the boundary layer theory and the separation condition uniquely 
specifies the strength and position of the nascent vortices and allows the 
separation points to move along the cylinder as required by the interaction 
between the wake and .the boundary layer. It is of importance to note that 
the no-slip condition is satisfied only at discrete time steps. Efforts 
to satisfy the said condition in a more continuous manner are counterproduc- 
tive because the decrease in At calls for unacceptably large computation 
times. One of the essential features of the model is the bracketing of a 
suitable interval for At, which provides an adequate scale of discretiza- 
tion without excessive computation time. 

The method described above properly accounts for the vorticity generated 
on the forebody of the cylinder. However, the flow in the wake of the 
cylinder develops a boundary layer on the afterbody which generates 



A lengthy discussion may be made regarding the difference between the 
definitions of steady and time dependent separation criteria. However, a 
numerical analysis such as the DVM does not warrant the consideration of 
the distinction between the results stemming from the use of different 
separation criteria. 
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countervorticity. Procedures by which this relatively small amount of 
vorticity generation may be accounted for are far from simple. The velo- 
city distribution on the afterbody has strong spatial and temporal gra- 
dients due to the immediate proximity of the fluctuating wake. These 
gradients render useless the separation prediction methods described 
ear I i er. 

Davis C623 attempted to account for countervort ic i ty effects by using 
four separation points (two for the forebody and two for the afterbody 
boundary layers) but his method required the assumption of fixed separation 
points. Deffenbaugh and Marshall C47D approximated the rear shear layer 
effect by assuming that the boundary layers remain attached over the same 
fraction of a region of adverse pressure gradient as that of the forward 
boundary layers. However, in view of the complex velocity distribution 
on the afterbody, it is incorrect to apply the same boundary- I ayer con- 
siderations as were used on the forebody. It is safe to say that the 
numerical perturbations introduced into the model in attempting to account 
for the secondary effects of countervorticity may cause more error than 
would the neglect of the effect altogether. This consideration plus the 
Introduction of the circulation reduction mechanism (yet to be described) 
obviate the need for the separate accounting of the effect of the rear 
shear layers. 

E. CONVECTION AND REDISCRETIZATION MECHANISM 

The vortices emanating from a separation point may be convected in a 
number of ways. Three most promising schemes are those given by the 
following expressions: 



zCt+At) = z(t-At) + 2q(t) At 


(29a) 


z(t+At) = z(t) + ^3q(t)-q(t-At)D At 


(29b) 


z(ttAt) = z(t) + q(t) At 


(29c) 
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where q(t) is the complex velocity of the vortex evaluated at the positions 
and times indicated. These schemes yield more or less the same accuracy 
and the use of one in lieu of the others is primarily governed by the sim- 
plicity of the method and computation time considerations. In the present 

t. 

analysis the simple Eulerian scheme given by Eq. (29c) is used because of 
Its simplicity and because the use of relatively small time intervals (At = 
0.125) make the method sufficiently accurate. Clearly, the results obtained 
with the three schemes cited above, or with others, cannot be compared 
without taking into consideration the effect of the time intervals used 
and the discretization process (if one is employed). All such convection 
schemes must assure, within the overall limitation of the model, that the 
vortex paths in regions of very small radius of curvature do not signifi- 
cantly deviate from those calculated with much smaller time intervals. 

It has been explained in connection with the discussion of the histori- 
cal development of the DVM that a I I investigators, with the exception of 
Fink and Soh C46H, let the vortices move in the form of vortex clouds and 
took special measures to eliminate the difficulties arising from vortex- 
vortex and vortex-cylinder proximity effects. Not to be entirely critical, 
the use of vortex clouds can yield force-transfer coefficients nearly 
identical to those obtained with the use of rediscretized sheets. However, 
one drawback of the vortex clouds is that it does not allow the formation 
of spiralling vortex sheets as often observed experimentally. There are 
other features of the rediscretization which will become clear later. 

The vortex clouds stem partly from the accumulation of the logarithmic 
errors and do not in any way represent a turbulent vortex. In the early 
stages of the present investigation the vortex clouds were used to obtain 
some idea about the gross features of the asymmetric wake. Figure 10 is 
a sample computer plot of vortex clusters. It is not possible to identify 
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Fig. 10. Flow configuration without rediscretization. 




Fig. II. Flow configuration with rediscretization. 
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a spiralling sheet in Fig. 10, but the shedding of the vortices and the 
thinning of the connecting vortex sheets between two adjacent vortex 
clusters are clearly observable. 

The fact that the vortices must represent each segment of the shear 
layer as accurately as possible and thus avoid the buildup of logarithmic 
errors in their convection gave rise to the method of rediscretization or 
redistribution. A sample flow pattern obtained with rediscretization is 
shown in Fig. II where the point vortices are connected by line segments 
to denote the position of the sheet. The mathematical foundations of the 
concept have already been introduced. In the following, the details of 
the mechanism where the rediscretization is affected are described. 

Consider a segment of the vortex sheet at an instant t and assume that 
the vortices are connected by straight line segments (see Fig. 12). The 
distance between any two vortices is given by 




Fig. 12. Vortex sheet before rediscretization. 



6s = Iz - z ,1 (30) 

n ' n n- 1 ' 

The length of the sheet from the position Z| (core of the spiral) to the 
n-th vortex is given by 
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n 



S = y 6s. = y |z. - 2. , 

n Z-/ j Z-^ ' j j-l 

j=2 j=2 



(31 ) 



The total length of the sheet for N number of vortices becomes 



= Z -J 



(32) 



j=2 



The vorticity per unit length of the sheet at each is given by 



Y, = r,/Ss2 



The circulation per unit length, f^iay be tabulated as a function of s^ 
through the use of Eqs. (31) and (33). This tabulation is then used as a 
piecewise continuous function denoted by y = y(s). The function y(s) uses 
a linear interpolation of the tabulated values, i.e. 



^n'Vl 



y(s) = y , +(s-s , ) 

' 'n-l n-l s -s 



I 

n n- 1 



s , < s < s 
n-l — — n 



(34) 



y(s) represents circulation per unit length and is an invariant at a given 
time of the rediscretization process. 

The next step in the analysis Is the replacement of the existing un- 
evenly spaced N vortices with N evenly spaced new vortices. In general N 

ys 

may be different from N. Replacement of N vortices is accomp I i shed by 
starting at the core of the spiral, Zj, and placing vortices at intervals 
of 6s measured along the straightline segments connecting the original 
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vortex positions (see Fig. 13). It shouid be noted that when the method is 
appiied to the exampie depicted in Fig. i3, and are positioned such 
that I Z 2 ” I = "Ss and |z^ - Z 2 I + |z 2 - Z 2 I = 6s. In other words, dis- 
tances are measured along the straightline segments of the existing sheet. 




Fig. 13. Vortex sheet after rediscretization. 



The strengths of the new vortices are calculated through the use of the 

distribution function y(s) such that the circulation per unit length at a 
/\ 

new position, z^, is the same as for the original sheet at the said posi- 
tion, i.e.. 



I = y(o) 6S2 



r = y(s*) 



• 6s . i+6s , 
n+l n-l 






Y(Sjj) 6s 



I < n < N 



(35) 



where 



s* = (n- I ) 6s 



The final step in the basic procedure is to insure that there has been no 



net gain or loss of circulation during the rediscretization process. The 



circulation deficit is computed as 



Ar = E 




(36) 



and an equal portion of the difference is added algebraically to each vor- 
tex, i.e.. 







= r + 

n 



Ar 



(37) 



Two parameters yet to be specified are 6s and N. The present model 
does not change the number of vortices during rediscretization and thus 

/N ^ /S 

N = N. 6s and N are related by 



6s = S|^/(N - I ) (38) 



It is important to note that N is the total number of vortices in a detached 
sheet (The definition of a detached sheet is introduced later.) and the 
rediscretization process does not change the position of either the first 
or the N-th vortex. For the vortex sheets attached to a nascent vortex, 
however, N represents the number of vortices excluding the nascent vortex. 
The reason for this is that the conditions (strength and location of the 
nascent vortex) prevailing at or near the separation point should not be 
altered by any process once the nascent vortex has been placed in the flow 
in accordance with the no-slip condition. Thus, for an attached sheet the 
rediscretization process stops at the vortex immediately adjacent to the 
nascent vortex. 

It is appropriate at this time to give some idea about the actual 
magnitude of some of the parameters mentioned above: 
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(i) the strength of a nascent vortex = 0.05 to 0.15 

(il) N = 3 to 80 depending on the length of the sheet 

C i i i ) 6s = 0. I 

The above numbers are representative of those encountered for a stationary 
cylinder with a time step. At = 0.125. 

F. VORTEX SHEDDING AND CIRCULATION REDUCTION 
I . Vortex Shedding 

One of the least understood of the wake formation processes is vor- 
tex shedding where the sheet connecting the spiralling vortex to the sepa- 
ration point is cut. That portion of the sheet which remains attached to 
the cylinder begins to roll up into a new vortex spiral while the detached 
portion moves downstream to become part of the vortex street. The shedding 
process is accompanied by a number of Interactive phenomena such as the 
motion of the shed vortices, adjustment of the longitudinal and lateral 
spacing of the vortex cores, and the reduction of circulation through the 
exchange of oppositely signed vorticity. These phenomena are so inter- 
locked that attempts to explain the reasons for the occurrence of one be- 
come in turn questions to be explained. Evidently one possible way to 
resolve this dilemma is to have an exact solution of the problem. Numeri- 
cal models such as the finite difference technique, having failed to 
elucidate the interactive mechanisms, left researchers with the need to 
develop techniques which could be used in more phenomenological models 
such as the DVM. It is absolutely essential that one begin the analysis 
with a model and then examine its validity through its predictions. 

Inspired by flow visualization and wake measurements, Gerrard C63I! 
proposed that the growing vortex spiral draws the opposite shear layer 
across the wake when it becomes sufficiently strong. The approach of 
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oppositely signed vorticity cuts off further supply of circulation to the 
vortex. The approximate time of the cutting of the sheet is then considered 
as the beginning of the shedding process. This mechanism is sketched in 
Fig. 14. 

The wake comprised of such shed vortices continuously adjusts it- 
self in both lateral and transverse spacing. In general it is agreed that 
the wake may be divided into three parts. The first is the so called for- 
mation region - 5); the second, the stable region where the vortices 

exhibit the character isttcs of a fairly uniform vortex street (5 < Z < 

12); and the third, the unstable region extending beyond x larger than 
about 12 C64j. 

The lateral spacing of the vortex cores gradually increases within 
the stable region and, like the longitudinal spacing, depends on the 
Reynolds number. 

Even though considerable experiments have been conducted (particu- 
larly at low Reynolds numbers) to delineate the basic features of the wake, 
very few attempts have been made to measure directly the strength of the 
vortices as a function of x at a given Reynolds number. Admittedly, this 
is an extremely difficult undertaking. Several investigators (Schaefer 
and Eskinazi C64H, Griffin C65l], and Davies C66j) measured the longitudinal 
and transverse spacing of the vortices in the stable region and calculated 
indirectly the strength of the vortices (assumed to be identical) through 
the use of von Karman’s vortex street model. Consequently, the vortex 
strengths are approximate and highly dependent on the initial inputs. The 
most important conclusion resulting from such investigations is that the 
strength of the vortices decreases with increasing distance downstream. 

The question concerning how much of the circulation generated remains in 
a vortex has been the subject of discussion since the time of PrandtI. As 
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Fig. 14 Initiation of vortex-shedding sequence 



noted earlier Prandtl suggested that the net circulation remaining in a 
vortex is about 50 percent of that generated. Discrete vortex models used 
by Clements and Sarpkaya C68] give a reduction of about 15 percent. 

This value is signi f leant ly smaller than the estimates of 40 percent re- 
duction (Page and Johansen Cl03), 57 percent (Roshko Cl G) and 70 percent 
(Bloor and Gerrard C69!]). It should not be inferred from the foregoing 
that the vortices once having acquired a net circulation of about 40 per- 
cent continue to retain that circulation forever. The understanding of 
the physics of circulation reduction is of major importance in the devel- 
opment of any phenomonolog lea I model and is discussed in greater detail in 
the fol lowing section. 

2. Circulation Reduction 

All investigators using the DVM have recognized that greater atten- 
tion must be paid to the circulation reduction if the DVM is to produce 
results comparable to those obtained experimentally. In the early stages 
of the evolution of the present model the need for a circulation reduction 
mechanism was quickly reaffirmed from the magnitudes of the lift and drag 
forces. 

Before undertaking a detailed discussion of possible circulation 
reduction mechanisms. It Is advantageous to discuss the results obtained 
through the measurement of the vortex strengths. Schmidt and Tilmann [ 10 ^, 
through the use of an ingenious ultrasonic sound beam technique, measured 
the variations of the circulations of vortices without inserting a probe. 
Their measurements were made at various Reynolds numbers by changing both 
the velocity of flow and the diameter of the cylinder. Figure 15 shows T 
as a function of x and is representati ve of the results obtained by Schmidt 
and Tilmann. In their experiments the Strouhal number was about 0.21 and 
the mean velocity at the point of separation may be takerr as = 1.5 U« 
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Percent Circulation Reduction 




Fig. 15. Vortex strength versus distance [69, 70U. 




Fig. 16. Percent of circulation reduction versus distance [69, 70]. 
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Thus, the total circulation generated during the period T = l/f is approxi- 
mated by 



r = 




D 

U St 
o 



10.7 



In the observed formation region, I < x < 6, the vortices grow rapidly and 
at the same time lose some of their strength due to various circulation re- 
duction mechanisms. One can think of a net circulation reduction only for 
vortices which are in the stable and unstable regions. It must be noted 
that the ultrasonic technique used in Schmidt and Tilmann's experiments 
cannot yield accurate vortex strengths in the formation region because the 
path of the sound beam crosses over two oppositely signed vortices. In 
the stable region, however, only one vortex crosses the path of the sound 
beam at a given time. For the reasons cited above, part of the curve ex- 
tending from X = 10 to its experimentally reported maximum value of 40 is 
used to calculate the circulation reduction by dividing the circulation 
values by 10.7. The result is shown in Fig. 16. Also shown in Fig. 16 
are two reduction ratios computed similarly from Bloor and Gerrard's C69H 
reported vortex strengths at x values of 12 and 20. The circulation re- 
duction shown in Fig. 16 should be interpreted only as a reflection of 
the change of circulation with distance rather than as definitive values 
since such values depend on the Reynolds number, three-dimensionality ef- 
fects, and the methods of their evaluation. 

One or more physical mechanisms are needed to explain and to model 
the observed circulation reduction. Circulation may be reduced only by 
oppositely signed circulation. This hypothesis is somewhat mechanistic 
but is easily adaptable to numerical modelling. There is no question that 
in reality the entire mechanism may be related to the action of viscosity. 
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Leaving aside momentarily the complex relationship between the action of 
viscosity and circulation reduction, attention is concentrated on three 
sources of circulation exchange. It has been noted by the users of the 
DVM that the circulation may be reduced by the oppositely-signed circula- 
tion generated on the afterbody (countervort i ci ty ) . In accordance with 
this assumption, the discrete vortices which come closer than a preassigned 
distance to the cylinder and the real vortices which cross the cylinder 
boundary are annihilated. Experiments by Page and Johansen ClOD have shown 
that only about a 10 percent reduction is brought about by the entrainment 
of fluid bearing circulation of opposite sign from the flow aft of the 
separation points. Numerical experiments to be discussed in another sec- 
tion of this work, as well as those done by others, show that the circula- 
tion generated on the downstream face of the cylinder can account for no 
more than about l/5-th of the total reduction in circulation. This fact 
calls for other mechanisms which account for the remainder of the circula- 
tion loss. The second mechanism postulated for the exchange of oppositely- 
signed circulation is an integral part of the shedding mechanism. As pro- 
posed by Gerrard C63j, when the growing vortex acquires a sufficiently 
large strength. It draws the feeding shear layer of the other growing vor- 
tex across the wake, thereby inhibiting further growth of the vortex. 

The tail of the cut sheet is then absorbed by the still attached vortex on 
the opposite side of the wake. This mechanism, quite distinct from the 
effect of the boundary layer on the rear of the cyl inder, further reduces 
the strength of a given vortex. However, numerical experiments again show 
that the reduction due to this mechanism alone is far from sufficient to 
bring the circulation down to the observed or measured values. Consequent- 
ly, a third mechanism must be postulated with which not only the vortices 
in the formation region but also those further downstream experience 
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additional circulation reduction. The origins of the third reduction mech- 
anism may be partly in the ingestion or entrainment of fresh fluid into the 
shear layers, partly in the instabilities in the spiralling vortex sheets, 
partly in slow viscous reduction, and partly through some other as yet 
unknown mechanism. Regardless of the true origins of the reduction mech- 
anism, two facts remain undisputed: the circulation of a given vortex does 

decay with time and the vorticity diffuses through viscous and turbulent 
action; and, the DVM cannot predict even the most gross features of the 
real flows (for example, lift and drag coefficients and circulation distri- 
butions) without the use of a circulation reduction mechanism. 

Lastly, it is necessary to briefly discuss the effect of diffusion 
of vorticity on the apparent reduction of circulation. Diffusion means 
that the vorticity is spread over a larger area. This in turn means that 
larger numbers of suitably distributed point vortices for each segment of 
the vortex sheet rather than a single elemental vortex represent the corre- 
sponding distributed vorticity. Ordinarily, the point vortices (whether 
on a sheet or concentrated) are assumed to be placed at the center of 
gravity of the vorticity field. ^ This does not imply that the single vor- 
tex at the center of gravity wil I yield the same velocity at an arbitrary 
point on the cylinder or at the position of another point vortex as the 
distributed vorticity. As far as the cylinder is concerned, the said 
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effect becomes relatively small as the distance between the vortex and the 
cylinder increases. Consequently, the calculations are not noticeably 
affected by the diffusion of vorticity in the far field. Since diffusion 
is a matter of time, the vortices in the immediate vicinity of the cylinder 
do not undergo significant diffusion. The foregoing simply points out the 
difficulty of accounting for all direct or indirect decay mechanisms and 
the need for a suitable hypothesis which will account for the circulation 
decay and yet be compatible with the basic premises of the DVM. As is 
discussed in greater detail later, the present model accounts for the effect 
of rear shear layer, circulation cancellation brought about by the shedding 
mechanism, and finally by the cross-street circulation exchange. The last 
process is accomplished simply by transferring a certain percentage circula- 
tion between oppositely signed vortices at each time step. 
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IV. DETAILS OF THE NUMERICAL MODEL AND DISCUSSION OF RESULTS 



A. SPECIFIC DETAILS OF THE MODEL 

I . Typical Sequence of Calculations at an Arbitrary Time 

The most straightforward way to describe the particulars of the 
model is to examine the sequence of events that occur during a typical com- 
putational cycle. Figure 17 is a flowchart indicating the order in which 
the model executes the various phases of the analysis. Each of the events 
that occur during every computation cycle are described in the order of 
their occurrence. The asymmetry introduction and vortex shedding processes 
occur only under certain flow conditions and are addressed separately. 

At the beginning of each computation cycle the velocity distribu- 
tion on the forebody of the cylinder is calculated at one degree intervals 
from the most recent position of the forward stagnation point. The posi- 
tions of the upper and lower separation points are then determined through 
the use of Pohihausen's method as described earlier. This allows the 
strength and positions of the two nascent vortices to be determined accord- 
ing to the no-slip condition. 

The next step is the convection of each of the vortices by a simple 
Eulerian method where the velocity of each vortex is computed at its pres- 
ent location. The velocities of the nascent vortices are computed at the 
point of their introduction prior to their actual placement into the flow. 
This insures that a nascent vortex is convected with the velocity prevail- 
ing at the time and place of its creation. Subsequently, the lift and drag 
coefficients are determined through the use of Eq. (15). 
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Fig. 17. 



Flowchart of computation sequence. 
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The normal sequence of the calculations requires a test for the 
conditions to determine whether the vortex sheet shouid be cut and the 
attached vortex spiral be shed. Because of its importance and the need 
for special handling, it is discussed separately. 

In discussing the problems associated with the use of discrete vor- 
tices, it has been pointed out that the elemental vortices must be placed 
at the center of the sheet segment which they are supposed to represent. 

This led to the method of rediscretization. Thus, following the test for 
vortex shedding and if necessary, after the cutting of the vortex sheet, 
all attached and detached sheets are red iscret ized. In doing so, each sheet 
is separately rediscretized (see Fig. 18). The discretization of the de- 
tached sheets is handled in accordance with the steps described in Chapter 
II, Section E. The only exception to the said process is in the treatment 
of the attached sheets near the nascent vortex. It has been repeatedly 
pointed out that the conditions near the separation point must be disturbed 
as little as possible by red i scret i zat ion or any other flow of events. In 
accordance with this pr i nci p I e, the nascent vortex and its immediate 
neighbor are not allowed to change their strengths or positions. 

The next step in the sequence of calculations calls for circulation 
reduction and coalescence. It is simply hypothesized that every vortex in 
the flow field loses its strength in an amount proportional to its current 
strength. This hypothesis is based on heuristic reasoning, observations, 
and numerous calculations. Further justification of the hypothesis is made 
later. Here only the actual process incorporated into the computer code is 
descr i bed. 

During each computation cycle for t > 5, the strength of every vor- 
tex, except the nascent vortices, is reduced according to the following 
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detached sheet 
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Fig. 18. A sample wake configuration. 



expression 



r (new value) = (l-p/lOO) F (old value) (39) 

n n 

where p depends on the location = Re{z^} as depicted in Fig. 19. For 
t < 5, p is varied linearly from zero to its final value of unity. 




As discussed in connection with various circulation reduction mech- 
anisms, there are two other processes whereby the circulation is reduced. 
The first is the reduction resulting from the proximity of vortices to the 
cylinder. Whenever a point vortex comes closer to the cylinder than a 
radial distance smaller than Ar = 0.04, that point vortex is assumed to be 
ann i h i I ated and removed from the flow field. The part of the sheet con- 
nected to the nascent vortex is excluded from this process. The annihila- 
tion of circulation in this manner results only in a small reduction in the 
total circulation and cannot account for a 40 to 50 percent reduction in 
the strength of the vortices in the stable region. The previous applica- 
tions of the DVM which considered this annihilation mechanism together with 
that resulting from the occasional proximity of oppositely-signed vortices 
(rediscretization not used) obtained only a mere 15 percent reduction in 
circulation C67, 68j. This fact was the first indication of the need for 
an additional circulation reduction mechanism. 
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The second of the two mechanisms concerns the entrainment of the 
tail of the detached sheet into the vortices across the wake. For the sake 
of discussion, assume that an attached vortex sheet is severed at a suit- 
able point. The details of how the cutting is accomplished is discussed 
later. The tail of the detached vortex sheet is drawn into the region be- 
tween two attached sheets because of induction. If left alone that tail 
eventually develops complex kinks and spirals and makes the rediscretiza- 
tion process unnecessarily complex. As proposed by Gerrard C63], the cir- 
culation in the tail gradually ends up in the attached sheet across the 
wake. In the model a similar procedure is employed. That portion of the 
tail which penetrates into the mixing region beyond a line tangent to the 
two attached vortex spirals (upstream of line AA' shown in Fig. 20) is re- 
moved from the flow field, and its circulation is uniformly subtracted 
from that of the sheet across the wake. Calculations which keep track of 
this cross-wake circulation reduction show that approximately five percent 
of the circulation of the detached sheet is transferred across the wake 
by this mechanism. It should also be noted that the position of the line 
severing the tail of the detached sheet is not of particular significance. 
One may have used other suitable and defensible procedures to accomplish 
the same purpose. 

2. Initiation of Asymmetry and the Vortex Shedding Processes 

The first of these two events occurs only once and the second every 
time period a vortex is shed. As such, they constitute special events in 
the computational sequence and are discussed separately, 
a. Asymmetry Introduction 

In nature the inception of asymmetry and the shedding of vor- 
tices are intrinsic characteristics of flow about bluff bodies. The causes 
of asymmetry and the subsequent vortex shedding are closely linked to the 
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Circulation reduction associated with vortex shedding process. 



instabilities in the shear layers. These instabilities may be brought about 
by the slight nonuniformities in the ambient flow, randomly distributed tur- 
bulence, asymmetry of the body, etc. It appears that the symmetric state 
of separated flow is not its most stable state. 

Whether the flow analysis is by the DVM, finite difference, 
finite element or whatever, asymmetry must be artificially introduced. 
Otherwise, the flow field remains forever perfectly symmetric and the vortex 
shedding process never begins. Figure 21 illustrates the results of the 
present model when asymmetry is never forced to develop. 

Although the onset of asymmetry and the subsequent evolution 
of the wake have been observed and recorded by many researchers, little is 
understood about the controlling mechanism. The question of whether the 
asymmetry begins in the wake and is fed back to the boundary layer, or 
vice versa, is unknown. Perhaps both occur simultaneously. Regardless, 
the process is not well enough understood to suggest any preferable way 
to artificially initiate its occurrence. Consequently, the method utilized 
in the present model is a result of heuristic reasoning based on well 
known experimental observations and numerous numerical experiments with 
the model. Before describing the details and results of the various at- 
tempts, it is worthwhile at this point to relate some of the general ob- 
servations, or lessons learned, from the investigation. Firstly, whatever 
method of asymmetry introduction Is chosen, the steady state results of 
the model are identical; i.e., the model has no memory of the conditions 
which initiated the shedding of the first vortex. Secondly, as stated 
earlier, asymmetry need be introduced only for the shedding of the first 
vortex. Following the shedding of the first vortex, the shedding process 
continues under its own impetus. Thirdly, every asymmetry mechanism in- 
volves a trade-off between the intensity and the duration of the perturbation 
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Fig. 21. Evolution of flow without asymmetry introduction. 
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applied. Associated with the question of duration is the decision of when 
to initiate the mechanism. Fourthly, there is an optimum time during 
which the symmetric fiow is most susceptibie to the onset of asymmetry, 
if asymmetry is introduced during this period, the intensity and duration 
of its appiication can be minimized. Finaiiy, of the uni imited number of 
mechanisms that could be devised, it is best to choose the simpiest method 
that produces the desired asymmetry. 

The interactive nature of the vorticity generation mechanism 
and the vorticity distribution in the wake suggest that asymmetry might be 
introduced either into the generation mechanism or into the wake itseif. 
Methods appiied to the vorticity generation method invoive, for exampie, 
artificiaiiy increasing the computed vorticity f i ux at one separation point 
Csay by iO percent) and simu itaneous iy decreasing the fiux at the other 
separation point. However, sample calculations with this procedure show 
that the method is not effective. Furthermore, the large asymmetry causes 
the computed separation angles to undergo unrealistic excursions from their 
mean positions. 

Introduction of asymmetry into the wake can be accomplished by 
either artificially altering the strength and/or position of the point vor- 
tices on the existing sheets. Changing the strengths of the vortices In 
one sheet relative to the other results in a sluggish response of the wake 
similar to that described above. The general ineffectiveness of altering 
vortex strengths to produce the downstream convection of one sheet relative 
to the other suggests the more direct approach used in the present model. 
The rationale Is to initiate the downstream convection of one sheet by 
actually displacing it a small amount each computation cycle over a short 
period of time so that when the mechanism is discontinued the sheet con- 
tinues to move without external influence. 
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Numerous numerical experiments have shown that a gradual appli- 
cation and buildup of the displacement increment serves to provide a smooth 
transition from symmetric to asymmetric conditions. The initiation time 
and duration of application are based primarily on two considerations. 
Firstly, this time interval should correspond to the period during which 
the flow is most susceptible to asymmetry. Secondly, the initiation of the 
mechanism should begin after the initial peak in the drag curve. Since the 
cause of this drag overshoot is the symmetric growth of the vortices, the 
asymmetry mechanism should not be allowed to influence this symmetric de- 
velopment. 

In the following the specifics of the asymmetry introduction 
method are presented. During each computation cycle within the interval 
5.0 < t < 9.0 the vortices of the upper attached sheet are displaced down- 
stream an amount given by (see Fig. 22) 

Ax = 0.a2 {1.0 - cosCiT(t-5)/4D} (40) 




Fig. 22. Variation of the asymmetry perturbation Ax with time. 
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Several points are worthy of discussion. Firstly, the displacement given 
by Eq. (40) is not unique and any other form of displacement function could 
have accomplished the same purpose. The particular choice, depicted by 
Fig, 22, introduces the disturbances gradually and continuously without 
shocking the flow field. Lastly, the total displacement applied to the 
sheet during the time interval 5 ^ t £ 9 is only 8 percent of the radius 
of the cy I i nder. 

b. Vortex Shedding Process 

Vortex shedding is understood as the mechanism whereby the feed- 
ing layer is cut. The specific details of the conditions which dictate the 
separation of a spiralling sheet from its feeding layer are not known and 
have yet to be determined experimentally. The numerical model is not in 
need of such a cut were it not due to the I imitations of computation time. 

If the flow field is allowed to evolve asymmetrically, the flow pattern 
shown in Fig. 23 results. This figure shows not only the evolution of the 
doubly-connected vortex sheets but also the drawing- in of the tail of the 
shedding vortex into the mixing region. The continuation of the calculations 
in this manner are prohibitive even with the largest computers. However, 
it is instructive in two ways to examine the evolution of the wake in this 
manner. Firstly, the uncut sheet begins to change the sign of its curva- 
ture as the lift force goes through its extremal (see point A of Fig. 23). 

The change of curvature is in fact the beginning of the drawing-in process 
of the shear layer across the wake. Subsequently, a strong inflection and 
kink develop drawing the sheet closer into the formation region. Simul- 
taneously, the lift force increases towards its next extremal. The gradual 
nature of this process does not permit one to identify a specific moment 
at which the shedding should commence. Thus, the examination of other flow 
features exhibiting stronger reaction to the precise moment of shedding are 
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necessary. Numerous calculations with symmetric and asymmetric evolution 
of the wake and with and without the rediscretization of the shear layers 
have shown that the sheet to be cut bears its minimum circulation per unit 
length at a time immediately after dF/dt goes through a minimum (see Fig. 
24). In other words, the connecting sheet is weakest and, therefore, 
most susceptible to breaking when its circulation is minimum. It is on 
the basis of this fact that the sheet is cut immediately after the rate 
of circulation reaches its minimum. At the time of cut, the lift force 
is about 50 percent of its preceding extremal (see Fig. 25). 

The actual cutting process is accomplished by removing a single 
point vortex, for only one time step , from the sheet at a distance s = 0.4 
(Fig. 26). Subsequent calculations separately rediscretize the detached 
sheet and the part still connected to the nascent vortex. As the flow 
evolves, the attached sheet continues to grow rapidly to form a new spiral 
while the detached sheet moves further downstream. At each time step the 
computer code tests the condition as to whether the new connected sheet 
should be cut. At the time the said condition is satisfied, the vortex 
sheet is cut as described above. It Is clear that the shedding of the vor- 
tices and the cutting process alternate with a time period equal to half 
the Strouha I period. When a vortex sheet is cut, say at the lower side of 
the cylinder, the sheet which was previously detached from the top part of 
the cylinder is coalesced into a single vortex at its center of gravity. 

At the time of its coalescence, the core of the vortex spiral has moved a 
distance at least seven radii downstream. The calculations have conclu- 
sively shown that the process of coalescence at this distance from the 
cylinder does not alter the continuity of any variable (e.g., lift, drag, 
circulation, etc.). 
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dr/dt prior to and after the cutting of the sheet. 




Fig. 25. Variation of lift during a cutting cycle. 
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Fig. 26. Cutting of a vortex sheet. 
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B. DISCUSSION OF RESULTS AND SENSITIVITY ANALYSIS 

In the following, the results obtained with the numerical model (re- 
ferred to hereafter as the standard run) are presented. Following this, 
the variation of some of the important kinematic and dynamic parameters 
to the variation in others is discussed. The purpose of these discussions 
is to partly establish the stability of the results, partly to relax the 
range of variability of the nondisposable parameters, and partly to estab- 
lish the cause and effect relationships. At this time, no distinction 
is offered concerning the choice of parameters in the model and the calcu- 
lated results. Furthermore, no distinction is made between prediction and 
inspired hindsight. 

I . Results Obtained with the Standard Run 

The standard run is computed using At = 0.125, p = I .0 (see Fig. 19) 
and with the asymmetry introduction and the shedding mechanisms previously 
described. 

The evolution of the wake over one arbitrarily chosen cycle is shown 
in Figs. 27 and 28. Figure 27 depicts the near wake region and Fig. 28 
shows the far wake vortex street. 

The variations of the lift and drag coefficients with t are shown 
in Figs. 29 and 30. The drag coefficient reaches a maximum at t = 4.2 and 
then rapidly decreases to its ultimate mean value of about 1.2. The fluc- 
tuations of Cjj have a frequency twice the vortex shedding frequency and an 
amplitude of about 0.06. The overshoot of the drag coefficient is a direct 
consequence of the rapid accumulation of vorticity in the symmetrically 
growing vortex spirals. Following the shedding of the first vortex, 
drops sharply and the I i ft coef f i c lent increases, ultimately reaching a mean 
amplitude of about 0.65. The slight variations in the oscillations of 
and C|^ come from various sources. The vortices absorbed during the cross- 
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Fig. 27. Evolution of flow in the near wake. 
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Fig. 27 (con't). Evolution of flow in the near wake. 
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Evolution of flow in the far wake and the vortex street 
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Fig. 29. Drag coefficient versus time. 
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Lift coefficient versus time. 



wake cancellation mechanism may be one more or one less during two different 
cycles. The increase of the otherwl se finite length of the wake does not 
permit the vortex sheets to follow exactly Identical paths from one cycle 
to another. Regardless, the said perturbations In Cp and are quite 
small and well within the accuracy of the numerical model. 

The amplitude of the lift coefficient is certainly within the 
range of lift coefficients reported in the literature. However, it is not 
possible to find a truly two-dimensional laboratory experiment to compare 
the calculated results with greater definitiveness. More is said later 
about the dependence of the amplitude of Cj^ on other parameters, in par- 
ticular on the circulation reduction mechanism. 

Figure 31 shows a comparison of the calculated and measured drag 
coefficients. The results are strikingly similar and may be taken as a 
first indication of the workings of the model. 

The oscillations of the separation points in terms of the angles 
measured from 6=0 (not necessarily the instantaneous position of the stag- 
nation point) are shown in Fig. 32. The separation angle for a given sepa- 
ration point oscillates nearly sinusoidally with an amplitude of A6 = 3 
degrees about a mean value of 0^ = 77 degrees. The frequency of oscilla- 
tion of 6^ is identical to that of the shedding of vortices. 

The stagnation point does not remain at 6 = 0 even In steady flow 
about a stationary cylinder. Figure 33 shows the oscillation of the stagna- 
tion point. Similar oscillations have been reported by Dwyer and McCroskey 
C7G. The separation angles relative to the instantaneous position of the 
stagnation point are shown in Fig. 34. A comparison of Figs. 32 and 33 
show that the oscillations of the separation and stagnation points are 180 
degrees out of phase. In other words,. if the stagnation point moves up, 
say to 6 = 2, the separation point at the upper half of the cylinder moves 
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Comparison of predicted and measured drag coefficients versus time. 




101 



Fig. 32. Separation angle versus time. 
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Fig. 33. Stagnation angle versus time. 
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Fig. 34. Separation distance versus time. 



backwards from a mean position 0^ = 77 to 0^ = 74 degrees. This fact is 
also in conformity with previously reported results C7l3. 

The rate of change of circulation at both the upper and lower sepa- 
ration points is shown in Fig. 35. Because the generation and shedding of 
circulation are the essence of the entire bluff body problem, dP/dt has a 
special significance among all other results predicted by the numerical 
model. A cursory examination of this figure shows that there is a perfect 
similarity of the rates of circulation of the upper and lower nascent vor- 
tices. Secondly, the variations in dP/dt from one cycle to another as t 
increases from 40 to 200 are very minor, indicating that the vortex street 
may be regarded as infinitely long for all intents and purposes for t 
larger than about 40. 

The mean value of |dP/dt| is about 1.2. Noting that the mean base 
pressure is related to |dP/dt| by C723 

Cp^ = 1.0-2 |dr/dt| (41 ) 



The present calculations yielded Cp^^ = -1.35 which compares favorably with 
that reported by Roshko Cl G at high subcritical Reynolds numbers CCp|^ 
(experimental) = -1.36 for 20,000 < Re < 100,0003. 

The amplitude of oscillation of dP/dt is about 0.18 which indicates 
that the velocity at the separation point varies from 1.41 to 1.65. 

All of the parameters cited above, namely lift and drag coefficients, 
separation and stagnation angles, and dP/dt oscillate with a frequency 
identical to the vortex shedding frequency f. The Strouhal number, St = 
fyO/U is found to be St = 0.205. The experimental values ot the Strouhal 
number are reported to lie between 0.195 and 0.210 for the Reynolds number 
range from 2,000 to 400,000 (see Fig. 36). Thus, the predictions of the 
model are In conformity with the measurements in the range of Reynolds 
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Rate of circulation generation versus time. 
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Fig. 36. Strouhal number versus Reynolds number C73]. 
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numbers where the DVM is applicable. As is noted later, the Strouhal number 
is the least sensitive of all the parameters to the variations of the non- 
disposable parameters. 

One of the often discussed characteristics of the vortex street, 
particularly in the stable region, is the relative spacing of the vortices. 
Von Karman's classic stability analysis gave = 0.281 (see Fig. 37). 

Since von Karman, the said ratio has been determined experimentally by a 
great number of Investigators and it has been found that varies from 

about 0.19 to 0.3. 

The lateral and transverse vortex spacings obtained with the 
standard run are shown in Fig. 37. The asymptotic values of and h^ are 
well within the experimentally determined values C74l]. In the stable region 
of the street, the lateral spacing decreases somewhat and the transverse 
spacing gradually increases from about 0.63 to its asymptotic value of about 
1.05. The ratio in the asymptotic range is seen to be 0.23, a value 

which compares well with those reported experimentally C?4]. Evidently, 
it is not meaningful to think of vortex spacings within the formation 
region. 

All discrete vortex models of separated flow about bluff bodies 
indicate the need for greater loss of circulation in the wake in order to 
bring a closer correspondence between the calculated and measured separa- 
tion points, lift and drag coefficients, dP/dt, etc. For example, Deffen- 
baugh and Marshall C473 obtained a mean separation angle of 67 degrees and 
a mean drag coef f ic lent of 0. 9 even though the smaller separation angle 
should have yielded larger drag than that found experimentally due to the 
Increased wake size. In the present model all conceivable circulation 
reduction mechanisms have been taken into consideration. Their variations 
with X or t may not be unique or in conformity with those obtained 
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Fig. 37. Longitudinal and transverse vortex spacings. 



experimentally but their characteristics follow all vortex shedding mech- 
anisms so far advanced. Figure 38 shows the circulation retained by a 
vortex sheet from its inception to its subsequent coalescence and convec- 
tion into the far wake. Also shown in this figure are the two experimen- 
tal points reported by Bloor and Gerrard C693 and the mean curve drawn 
through the data reported by Schmidt and Ti Imann C703 (see Fig. 15). The 
comparison of a 1 1 the results is better than expected in view of the fact 
that the experimental determination of circulation through direct or 
indirect methods is an extremely difficult and approximate task. In par- 
ticular, the ultrasonic sound beam technique used by Schmidt and Ti Imann 
yields the vortex strengths only approximately in the stable and unstable 
region of the vortex street. There are no measurements of the vortex 
strengths in the formation region. The foregoing does not imply that the 
particular variation assigned to the circulation loss mechanism (see Fig. 

19) is unique. In fact, the only major nondisposable parameter in the 
entire numerical model is the form assigned to the circulation loss. The 
need for circulation loss has been amply demonstrated. The form of the 
particular relationship can only be adjusted with hindsight by comparing 
all of the predictions of the model with those obtained experimentally. 
Because of the fact that the differences in measured and calculated circu- 
lations are further downstream in the wake, the use of the experimentally 
suggested p variation does not in any way affect lift, drag, separation 
angles, dF/dt, and the Strouhal number. This will become more evident in 
the course of the discussion of the sensitivity of the model. Calculations 
not reported herein with more exotic forms of p have ascertained the 
validity of the insignificance of variation of p beyond the stable region. 

2. Sensitivity Analysis 

This section deals with the effect of the variation of one parameter 
abouts its value used in the standard run on the remaining major parameters. 
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Numerical results 




CD 

O) 



no 



Vortex strength versus distance C69, 70]. 



All numerical analyses must use a finite mesh size or time interval. 
The magnitude of these quantities might affect the results in various ways. 

In the DVM a large At might bring in curvature effects particularly in flow 
regions where the radius of curvature of the spiral is small. This effect 
is brought about by the convection of the vortices along straight lines co- 
incident with their instantaneous velocities. Thus At should be kept reason- 
ably small to allow the vortices to follow their natural paths. A very 
small At, on the other hand, requires prohibitively large computation times 
and defeats the purpose of the numerical experiments. 

The standard run was made with At = 0.125. Additional runs were 
made varying At from 0.1 to 0.2. The resulting drag and lift coefficients 
are presented in Figs. 39 and 40. The differences in the two coefficients 
for ali values are negligibly small and certainly within the range of 
accuracy expected from numerical calculations. Other parameters such as 
separation and stagnation angles, circulation loss, Strouhal number, etc. 
remain essentially identical to their corresponding values obtained with 
At = 0. i25. 

The role played by the loss of circulation in both iaboratory and 
numerical experiments has already been discussed. Here attention is devoted 
to the seiection of the particular value of p and the effect of its varia- 
tion on the remainder of the predictions. Several hundred hours of calcu- 
lations have shown that one can arrive at an approximate value of p by 
comparing the predicted and measured values of the vortex strengths and the 
lift and drag coefficients. None of these parameters, however, exhibit the 
same degree of variation with p at corresponding times and downstream dis- 
tances. For example, the comparison of the vortex strengths in the stable 
region, where most of the measurements are made, is not very meaningful 
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Fig. 39. Drag coefficient versus time for At=0. I , 0.125, 0.15 and 0.2. 




Fig. 40. Lift coefficient versus time for At=0.l, 0.125, 0.15 and 0.2. 
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since the effect of such vortices on the lift and drag coefficients is very 
small. Thus, it was decided to choose the p value such that the terminal 
value of the drag coefficient conforms to that measured in the subcritical 
region. No other experimental data was used in the analysis. Numerical 
experiments with various values of p have shown that p = 1.0 yields - 
1.2. Having determined a value of p, it is necessary to examine its varia- 
tion on a 1 1 of the predictions of the model. 

The drag coefficients obtained with p = 0, p = 0.5, p = I, and 
p = 2 are shown in Fig. 41. The insert in Fig. 41 shows the variation of 
Cp with p in the asymptotic range of Cj^. Evidently a 100 percent change in 
p about p = I caused only about 10 percent change in C^. Consequently, the 
tuning of the model through the use of = 1.2 is not as critical as it 
may have appeared at first sight. Figure 41 also shows that a circulation 
loss mechanism is necessary (see curve for p = 0) if the numerical model 
is to yield results in conformity with the experiments. As noted earlier, 
the need for circulation loss is also evidenced by the vortex-strength 
measurements (see Fig. 38). Thus, the important question in the evolution 
of the model is not whether there should be a circulation loss but rather 
what should the mechanism and magnitude of this loss be. 

The variation of the lift coefficient with p is shown in Fig. 42. 

It is immediately apparent that the lift coefficient is far more sensitive 
to p than Cj^. Regardless of its sensitivity, however, remains well 
within the experimentally reported values. One could fine-tune the value 
of p by matching the measured and calculated values of Unfortunately, 

thi.s is not possible since there are no reliable Cj^ values, particularly in 
the region of high subcritical Reynolds numbers. Consequently, no attempt 
is made to slightly increase or decrease p from p = I to obtain a specific 
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Fig. 41. Drag coefficient versus time for p=0.0, 0.5, 1.0 and 2.0. 
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Fig. 42. Lift coefficient versus time for p=0.0, 0.5, 1.0 and 2.0. 
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value. The importance of Fig. 42, rather, lies in the fact that the 
numerical model clearly shows the sensitivity of to circulation loss and 
offers an indirect explanation of the difficulties encountered in the 
laboratory experiments. The findings of the model can also be taken as a 
means of reducing the transverse force by artificially increasing the 
circulation dissipation in vortices (e.g., use of fairings, wire screens, 
etc.). 

The Strouhal number displays a relatively small variation with p 
Csee Fig. 43). This exp ia ins the good correspondence obtained between 
measured and caiculated values of St even when the circulation reduction 
does not exceed i5 percent (Clements C67H» Sarpkaya C6SD)* 




Fig. 43. Variaticn of the Strouha i number with p. 

it is of interest to point out the relationship between the strength 
and motion of vortices in the near wake and the sma i i variations in the 
Strouhal number. When the detached vortex is relatively stronger (p < I), 
then it iingers ionger near the cyiinder by virtue of the veiocities in- 
duced on it by the equa I iy.-strong vortex across the wake. This in turn in- 
creases the period of vortex shedding or decreases the Strouhal number. 

For weaker vortices (p > I) the inverse is true. 
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In the following, the sensitivity of the lift and drag coefficients 
and the Strouhal number to the variations in the amplitude of p has been 
discussed. In doing so, the shape of the dissipation function has not 
been changed. From a scientific point of view, the next to consider is 
the effect of the variation of the dissipation function on the major predic- 
tions of the model. For this purpose an extensive study was undertaken to 
acquire a deeper understanding of the dissipation particularly in the 
near wake region of the cylinder. In doing so, not only the terminal 
value of the drag coefficient but also other parameters such as the lift 
coefficient, separation points, and the amount of circulation retained by 
the concentrated vortices were carefully observed. Here only the high- 
lights of these observations can be described. The use of a dissipation 
function such as p = cr/x (for x > I) where a is a constant shifts the re- 
gion of larger dissipation closer to the rear of the cylinder. For example, 
the results obtained with o = 2.0 have shown that the lift and drag coeffi- 
cients show no more than a 10 percent variation from those obtained with 
the standard model. The separation points remained at almost exactly the 
same average values. The percent circulation retained by the concentrated 
vortices increased fpr obvious reasons. In fact, the relative vortex 
strengths at comparable disiances were somewhat larger than those measured 
by Gerrard C243. The Strouhal number decreased from 0.205 to 0.193, in- 
dicating that not only the vortices in the formation region but also those 
further downstream determine the transport velocity of the vortices and 
hence the Strouhal number. 

The experience gained from this exploratory investigation has shown 
that there should be a fairly large and uniform dissipation in the near 
wake, say within the formation region, and then decrease gradually with 
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di stance. However, wi thout. caref u I measurements of the vortex strengths 
as a function of time or distance it may never be possible to determine 
the precise shape of the dissipation function. 

To further determine the sensitivity of the model to dissipation and 
the relationship between someofthe most important parameters (e.g., 
separation angles, vortex strengths, and Strouha I number), ari additional 
series of numerical experiments have been performed. Specifically, the 
separation angle was determined not in the manner described in the stan- 
dard model but rather by suitably extrapolating the velocity profiles 
beyond the point of maximum velocity. The use of such an extrapolation 
alleviated the influence of the nascent-vortex proximity effect on the 
separation point and resulted in larger separation angles (0^ = 81 degrees) 
When the shape of the dissipation function was kept exactly as in the 
standard run, the drag coefficient reached an ultimate value of about 1.3 
and the lift coefficient a value of about 0.8. The Strouhal number de- 
creased to about 0.19. It was thus apparent that the value of p in the 
standard dissipation should be somewhat increased, say to p = 1.3. Such 
an increase in p would decrease both and and increase the Strouhal 
number without affecting the separation and stagnation points. This 
numerical experiment suggests that it is possible to match the observed 
and calculated separation points and resistance to fluid motion with one 
nondisposable parameter. No attempt has been made to find a precise value 
or shape for the dissipation function because of the primary concern with 
the interrelationships between various parameters. The fact has been 
demonstrated that the physical mechanism governing the flow is not affected 
by the actual value of p or and C^. It Is evident from the foregoing 
that the numerical modelling of a turbulent wake does not require a 
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I 






dissipation based on turbulent diffusion and that diffusion could be 
quantified only by additional measurements which will complement the 
numerical experiments. As is often the case in many aspects of fluid 
mechanics, the question seems to have been reduced once again to the 
understanding of the physics of turbulence. 
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V. SUMMARY OF THE RESULTS 



The numerical model based on the rediscretization of vortex sheets and 
the mutual interaction between the boundary layers and the wake yielded 
numerous results which may be summarized as follows: 

1. A discrete vortex model of flow past a stationary cylinder may be for- 
mulated with only one nondisposable parameter. Experiments and the numeri- 
cal results show the necessity of circulation reduction but not necessarily 
the exact mechanism by which it is accomplished. 

2. Based on a postulated circulation reduction mechanism, all predictions 
of the model are in conformity with those obtained experimentally. The 
fact that all and not only a few of the parameters agree with the experi- 
mental measurements adds further credibility to the assumptions made and 
to the power of prediction of the model. 

3. A concerted effort is made to establish a balance between simplicity 
and accuracy. To this end, simple linear relationships have been preferred 
in lieu of more complex ones. 

4. One of the most important aspects of the numerical model is its flexi- 
bility. Practically all parameters have a wide stable range and allow 
the user to perform numerical experiments to seek a physical understanding 
of the intricate relationships between two or more parameters. 

5. The demonstration of the applicability of the discrete vortex model 
to the analysis of flow about a bluff body with mobile separation points 
shows that similar models can be constructed with ease for bodies with 
fixed separation points. This allows one to study such problems as vortex 
shedding from ship bilge keels, stall over aircraft, fuel slushing in large 
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containers, etc., in naval hydrodynamics and aeronautics. 

6. The model, as it stands, can be used in the analysis of a great many 
f low-structure interaction problems. Of particular interest is the analysis 
of a circular cylinder undergoing forced or self-excited transverse oscil- 
lations. Clearly, the numerical experiments with oscillating cylinders 
can not only shed additional light on the mechanisms which bring about such 
oscillations but also increase the universality of the model. 
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VI. APPLICATION OF THE DISCRETE VORTEX MODEL 
TO SELF-EXCITED TRANSVERSE OSCILLATIONS 



A. INTRODUCTION 

Numerous experiments have shown that (see Parkinson C75] for an ex- 
tensive review) when the natural frequency of a bluff body is close to 
the vortex-shedding frequency, the body may undergo self-excited trans- 
verse oscillations. The primary consequences of these oscillations are 
as follows: (i) The vortex shedding frequency locks on to the natural 

frequency of the body, (ii) The amplitudes of oscillation and the 
transverse force increase rapidly up to the point of perfect synchroni- 
zation where the maximum amplitude is attained. Subsequently, the oscil- 
lations decrease either rapidly or abruptly. The abrupt drop in amplitude 
is a consequence of hysteresis. The reasons leading to the hysteresis 
are not yet known, (iii) The transverse force leads the displacement 
by a phase angle which increases gradual ly at first to about 50 degrees 
and then rapidly to about 130 degrees, (iv) At the end of the synchroni- 
zation region, the vortex shedding frequency jumps to that governed by 
the Strouhal relationship. 

The studies concerned with the understanding and prediction of the 
synchronization phenomenon may be classified in three categories. The 
first and most extensive category has been the experimental observations 
and measurements with forced or self-excited oscillations of circular 
cylinders, cables, and a small group of other bluff bodies. The second 
category of investigations dealt primarily with the analysis of the 
structural aspects of the equations governing the vibrations by making 
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ad hoc assumptions about driving fluid forces. Finally, the third cate- 
gory concentrated on the understanding of the flow field and the nature 
of the fluid forces. It has not yet been possible to understand the 
coupled structural and fluid mechanical mechanisms to explain the reasons 
leading to the self-excited oscillations. 

Mathematical models based on the assumption of a nonlinear oscillator 
have been developed. Among those, the most noteworthy is that proposed 
by Hartlen and Currie C76]. Their model predicts some of the observed 
features of the vortex-induced oscillations. However, the model requires 
a number of nondisposable parameters and gives no clues about the cause 
and effect relationships. It has been realized, as suggested by 
Parkinson, C75j among others, that if any understanding of the phenomenon 
is to be achieved, theoretical developments and numerical experiments 
must be compared with laboratory experiments to guide and complement each 
other. It is with this suggestion in mind that the discrete vortex model 
has been applied to the prediction of the characteristics of self-excited 
osci I lat ions of a circular cylinder. In the following, first the 
description of the mathematical foundations of the problem, then the 
specific details of the numerical procedure, and finally the discussion 
of the results and a summary are presented. 

B. FORMULATION OF THE PROBLEM 

The equation of motion for an elastically-mounted and linearly-damped 
cylinder of unit length (see Fig. 44) may be written in dimensional form 
as 

• 2 2 

my t 2m 0)^ Cy + y = U c (42) 
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where m represents the mass of the cylinder; y, the displacement; oj^, the 
natural circular frequency; the damping ratio; p^, the fluid density; 
c, the radius of the cylinder; and C^, the instantaneous value of the 
lift coefficient. In Eq. (42), y and y represent respectively the first 
and the second derivatives of y with respect to time. 




Fig. 44. Elastically-mounted, linearly-damped cylinder. 



Defining f as the vortex shedding frequency for the stationary 
o 

cylinder, as the corresponding Strouhal number, and 
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Eq, (42) may be written in normalized form as 



2ircS ITS - , » 

fi + n + ( — -)^ n = 2tt^s ^a c, 



(44) 
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Because of the need to distinguish the stationary and oscillating vortex 
shedding frequencies, and hence the Strouhal numbers, the previously and 
commonly used notation of St for the Strouhal number is modified as 
denoted above. 

The evaluation of Eq. (44) requires the instantaneous value of the 
lift coefficient. In the present model it is evaluated through the use 
of the complex velocity potential modified to take into account the motion 
of the cylinder. Thus, one has 
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in which denotes the instantaneous position of the center of the cylin- 
der and z^, its velocity (see Fig. 45). The use of Eq. (45) together with 
the generalized Blasius theorem yields 
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Fig. 45. Coordinate system for oscillating cylinder. 
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A comparison of Eqs. (46) and (15) shows that is replaced by the rela- 



whlch account for the forces arising from the velocity and acceleration 
of the cy I i nder. 

The calculations for the boundary layer and the separation points re- 
quire the evaluation of the tangential velocities relative to the cylinder. 
At an arbitrary time t and angle 9 (measured as shown in Fig. 46), the 
tangential velocity is given by 



tive position in the latter equation and two new terms are added 



10 • 

u^ = Re {Cq(z^ + e ) - z^3 e 



-i(e-ir/2). 



(47) 




o 



X 

Fig. 46. Relative tangential velocity on boundary of 
moving cy I i nder. 



As before, the velocity at any point in the flow field is given by 
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C. NUMERICAL PROCEDURE 



It is apparent from the formulation of the problem that the solution 
of Eq. (44) requires two parameters (? and a^) characterizing the mechani- 
cal system and the instantaneous fluid force imparted to the cylinder. 

The parameter a^, known as the mass parameter, may be readily calculated 
for a given cylinder of mass m, radius c, and fluid density p^ Usee Eq. 
C43)3* The damping ratio, ?, is a measure of the relative energy loss 
within the material due to molecular motion and does not include the so 
called fluid damping. All fluid forces acting on the cylinder are included 
in the proper specification of C|^ which represents the integrated effect 
of the surface-pressure loading by the fluid on the vibrating cylinder. 

The discrete vortex model developed for the stationary cylinder is 
used together with Eqs. (46) and (48) to determine the fluid forces acting 
on the cylinder. Equation (44) is integrated through the use of a fourth 
order Runge-Kutta method. Otherwise, no additional modifications are 
introduced into the procedures described In the previous chapter. In fact, 
the only specification made was the time at which the cylinder was allowed 
to respond to the fluid forces. Such a decision was necessary to avoid 
the coupling of two transient states resulting from the impulsive start 
of the flow and the transient response of the cylinder. This would have 
made the calculations unnecessarily complex. Instead, it was decided 
first to establish a flow with a sufficient number of vortices in the 
wake (e.g., by regarding the cylinder as stationary for t < 10) and there- 
by studying the response of the cylinder through a more or less steady 
uniform flow. The approximations signified by "more or less" are due to 
the fact in all numerical ca leu I at ions the vortex street is finite. In the 
present calculations, the cylinder is allowed to respond after t > 10, 
which was sufficiently large for all intents and purposes. 
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Mathematically speaking, there are no particular restrictions on the 

selection of the numerical values of 5 and a other than the fact that 

o 

they should be physically realizable. From the point of view of actual 
calculation, however, one must place some restrictions on ^ and a^ to 
limit the computation time. The smaller the value of ?, the larger the 
amplitude of oscillation. The larger the mass of the cylinder or the 
smaller the parameter a^ (keeping everything else constant), the larger 
is the transient time during which the cylinder oscillations asymptoti- 
cally approach their final vajues. Consequently, a very small value of 
C coupled with a rejatively small value of a^ would require many cycles 
of oscillation before the cylinder reaches its final amplitude. This 
makes the computation prohibitively expensive and time consuming. The 
foregoing arguments constituted one of the primary conditions in the selec- 
tion of ? and a^. The second and probably more important reason from a 
scientific point of view was the focusing of attention on the understand- 
ing of the physics of the phenomenon rather than on the solution of 
specific problems. With the foregoing considerations in mind, the ^ and 

a values shown in Table II were used in the calculations. The role of 
o 

the new parameter Sg = ?/a^ appearing in the table is discussed later. 

TABLE I I 



SET 




a 

0 




1 


0.02 


0.0125 


1 .6 


1 1 


0.04 


0.0250 


1.6 


1 1 1 


0.08 


0.0500 


1.6 


IV 


0.04 


0.0125 


3.2 
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D. DISCUSSION OF RESULTS 



I . I ntroduct ion 

The major part of the discussion of the results deals with the 
phenomenological predictions of the model and their comparison with those 
observed or measured experimental ly. For this purpose only one set of 
damping ratio and mass parameter is considered. The variations of the 
lift coefficient, relative amplitude, phase angle between the lift force 
and the cylinder response, and finally the ratio of the cylinder-oscilla- 
tion frequency to natural frequency is presented as a function of U^ 
defined by 



U 
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(49) 



In dimensional terms, U^ = U/a)^0 . The reasons for the choice of this 
particular parameter are twofold. Firstly, it is commonly referred to 
in the previous works C75j. Secondly, it is a measure of the velocity 
of flow past a cylinder of given diameter and natural frequency. 

Figure (47) shows the results obtained by Feng C77H with a lightly 

damped circular cylinder. As noted earlier, when the vortex shedding 

frequency, f , for a stationary cylinder approaches the natural frequency 

CU^ - 0.8 in Fig. (47)J, the vortex shedding frequency becomes nearly 

identical to the natural frequency of the cylinder. Furthermore, the 

frequency of oscillation of the cylinder, co^, nearly coincides with its 

natural frequency. As U^ increases towards 1.0, the amplitude of both 

the oscillation and the transverse force Increases rapidly while to/o) and 

n 

remain nearly constant at a value slightly under unity. As U^ In- 
creases further, both amplitudes drop either gradually or abruptly as 
illustrated in Fig. (47). The abrupt drop is a consequence of hysteresis 
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Fig. 47. Experimental results for self-excited transverse oscilla- 
tions of a cylinder C77]. 
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noted earlier. For U^. > l.l the frequency of vortex shedding, jumps 

to its stat ionary-cy I inder value, i.e., = u) . However, the cylinder 

o 

continues to oscillate at the natural frequency of the cylinder with very 

small amplitudes. One can, therefore, identify two major regions in the 

relationship between w , w , and u . For U smaller than about 0.8, to = 

V c n r V 

to and u - ui . In the range 0.8 < U < I . I , to < to and to = to = to . 

V cn ^ ~r~'vv vcn 

o o 

For l.l < U < 1.3, to = to and to == to . 

o 

The phase angle between the exciting force and the cylinder re- 
sponse undergoes dramatic changes in the region where 0.8 < U <1.1. 

The reasons leading to the amplification of the oscillations of the 
cylinder and the lift force, and the rapid increase of the phase angle 
from about 20 degrees to 100 degrees have not been understood because of 
the difficulty of measuring or observing the detailed kinematics of the 
flow field about the cylinder. Many experiments have essentially confirmed 
the general features of the results obtained by Feng C7?3. Most of the 
attention has been devoted to the understanding of the locking-on of the 
vortex shedding frequency to the natural frequency of the cylinder in the 
range of U values 0.8 < U <1.1. In conformity with the previous 
literature, this region is variously referred to as the synchronization 
region, locking-on region, capture region, or simply as the resonant region. 

The numerical results obtained with ^ = 0.02 and a^ = 0.0125 are 
shown in Fig. (48), again as a function of U^. It is immediately apparent 
that the results shown in this figure follow essentially the same trends 
as those in Fig. (47). Also shown in Fig. (48) are the corresponding 

values of co and ratios to /to* or U /U* where to* and U* are the values of 
o o o r r or 

to^ and when reaches its maximum value (point A). This point Is 
referred to as the point of perfect synchronization. 
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2. Mechanism of Synchronization 

Previous investigators C75U have uncovered a number of relatively 



unrelated facts regarding the kinematic and dynamic behavior of flow 
about a transversely-oscillating cylinder. The most notable of the ob- 
servations and measurements may be summarized as follows; 



(i) 


Oscillations increase the spanwise correlation; 


cm 


The strength of the vortices increases; 


(iii) 


The mean base pressure decreases and the mean drag coefficient 
i ncreases; 


(iv) 


Transverse force is amplified relative to its stationary-cylinder 
va 1 ue; 


Cv) 


The phase angle between lift and cylinder motion increases with 


(vi) 


The velocity range over which synchronization occurs increases 
with oscillation amplitude; 


(vi 1) 


Separation point excursion increases with oscillation; 


(vlii) 


Lateral spacing of the vortices decreases with increasing ampli- 
tude and is unaffected by changes in frequency of oscillation; 


(Ix) 


Wake-formation length decreases systematically with increasing 


(x) 


amplitude of vibration for U < U* and increases for U > U*; 

r r r r 

The wake contracts for U < U* and expands for U > U*. 

r r ^ r r 

The facts summarized above have been distilled from a number of 



investigations, not all of which have been conducted at identical Reynolds 
numbers or under similar experimental conditions. Consequently, there 
are certain features of the flow which have been observed to increase in 
one investigation and decrease in another. Whenever such conflicts arose, 
heavier weight was subjectively attached to the observation more frequently 
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reported. Perhaps the most interesting feature of all these isolated 
facts is that it has not yet been possible to establish the cause and 
effect relationship among them through a phenomonologica I ly and mathe- 
matically defensible mechanism. In the following, the mechanism suggested 
by the numerical experiments is described in detail. 

Consider an elastically mounted circular cylinder which is at its 
lowest displacement during a given cycle. Also assume that the stationary 
cylinder value of the vortex shedding frequency is slightly greater than 
the natural frequency of the cylinder. As soon as the cylinder begins 
to move upwards, the lower attached-vortex sheet lags behind the cylinder 
because of its inertia. In other words, the wake axis passing through 
the instantaneous center of the cylinder rotates clockwise (see Fig. (49) 
for t = 189). In the meantime, the velocity relative to the cylinder 
rptates clockwise an angle dictated by the velocity of the uniform flow 
and the instantaneous velocity of the cylinder. The rotation of the 
velocity vector coupled with the rotation of the wake axis moves the 
upper separation point further downstream and the lower separation point 
further upstream relative to that which would have existed on a stationary 
cylinder. The rate of circulation at the upper separation point increases 
rapidly and decreases at the lower separation point. This gives rise to 
a larger circulation per unit time than that for the stationary cylinder 
for the vortex towards which the cylinder is moving. Furthermore, the 
vortex spiral which Is attached to the upper separation point does not 
move out of the cylinder's way because of its inertia and consequently 
comes closer to the approaching cylinder. These two facts, namely the 
increase of circulation and the decrease of the relative distance between 
the cylinder and the upper vortex spiral, confirm the observation made by 
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Fig. 49. Wake configuration for self-excited transverse oscillations. 
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Fig. 49 (con’t.). Wake configuration for self-excited transverse oscilla- 
tions. 
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Davies C78] that "The growing vortex on vibrating cylinders seems to roll 

up more quickly and with a strength about 35 percent larger than that for 

a stationary cylinder." The increase of the vortex strength predicted 

by the current model will be discussed later. The phenomenon is, however, 

relatively more involved than the simplified introductory remarks would 

suggest. There is a third factor which enhances the strengths of the 

attached vortex sheets. This is directly related to the difference between 

the period of oscillation of the cylinder, T = l/f = l/f , and the time 

c c n 

required for the shedding of a vortex from a stationary cylinder, i.e., 

T^ - T^ . Each growing vortex about an oscillating cylinder continues to 
o 

be fed circulation not only at a higher rate but also for a longer time 
period. This further reinforces the strength of each growing vortex. 

It is now necessary to establish the link between the enhanced 
vortex strengths and the locking-on of the vortex shedding frequency to 
the frequency of cylinder oscillation. The results of the stationary 
cylinder have shown that the larger the strength of the vortices in the 
near wake the smaller their velocity relative to the cylinder. That is, 
strong vortices with their equally strong images tend to conserve their 
moment of circulation. Consequently, the vortices are shed at larger 
time intervals. This, in turn, leads to smaller Strouha I numbers. The 
decrease of S, or the increase of T ., brings T closer to T . This ten- 

V V ^ 

dency is reinforced every succeeding cycle until = T^ = T^ or until 
the effects which cause the changes in amplitude, period, lift force, 
and the phase angle disappear. At such time, a stable periodic state is 
reached. 

The foregoing explanation dealt only with the kinematic aspects 
of the fluid motion. In the final analysis, it is the balance between 
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the work done by the fluid on the cylinder and the work done by the dissi- 
pative forces (Internal friction) which determine the stability of the 
oscillations. If the net work becomes negative, then the oscillations 
decay and vice versa. The parameter which establishes the relationship 
between the friction mechanism and the given state of oscillation is the 
phase angle between the exciting force and displacement. It is given by 
C79D 

2? (ii/bi 

tan <j) = ^ (50) 

I - (w/o) ) 
n 

As approaches unity from < I , <{) increases rapidly, the rate of 

increase being a strong function of For a cylinder with very small 
Internal damping, (p increases from say 20 degrees to <|) = 150 degrees for 
a very small change in co/o)^. 

The relationship between the dynamics and the kinematics of the 
motion for the three states of oscillation and hence the reason for the 
sustained oscillations in the region of synchronization may be explained 
as follows. In the case of synchronization (see Fig. 50) the maximum 
velocity of the cylinder leads the maximum lift and the lift reaches zero 
at about the time the cylinder reaches its peak amplitude. In other 
words, the lift force opposes the cylinder motion for a short time period 
as the amplitude of the cylinder begins to decrease from its maximum 
displacement. If for any accidental reason the amplitude of the cylinder 
in one direction were to be increased while the cylinder was undergoing 
synchronized oscillations, the increased strength of the vortex growing 
ahead of the cylinder motion Increases the period, T^. This, in turn, 
decreases the phase angle between the force and the displacement which 
then i ncreases the phase angle between the force and the velocity. 
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Fig. 50. Lift coefficient and cylinder velocity versus time. 
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Consequently, the energy transferred from the fluid to the cylinder de- 
creases below that dissipated by the material damping. These sequences 
of events decrease the amplitude to the level where the energy extracted 
from the fluid over a cycle is just balanced by that dissipated by 
internal friction. 

The role played by all of the terms of the governing equations of 
motion and the energy transfer mechanism may be quantified as follows. 
Re-writing Eq. (44) 

271^$ , ttS « - ^ 

fi + — f) + ( ) q = 2ir S a C, (44 repeat) 

0 ) 0 ) o o L ^ 

o o 



and multiplying with q and integrating over a cycle, one has 



t+T 



t+T 



t+T 



t+T 



j qqdt + j q^dt ^’^dt = J C^q 

t t t t 



qdt (51) 



which reduces to 



t+T 



t+T 



2ttcS 



^^o f ’2 2 2 f 

— - I q dt = 2 tt S a / C, qdt 

o J o oj L' 



(52) 



In the case of a perfect balance between the energy input Cright hand side 
of Eq. (52)3 and the energy dissipated (left hand side), the above equality 
holds true. For an unstable system where the equality does not hold true. 



2 2 

E = 27T S a 
o o 



t+T t+T 

/ f " 



q^dt 



(53) 
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Thus, E = 0.0 corresponds to steady state with constant amplitude of os- 
cillation; E > 0.0 corresponds to ampi i f ication; and finally, E< 0.0 
corresponds to damping of oscillations. The variation of ri» h> H and C|^ 
in Eq. (44) as a function of time is shown in Figs. 5la-c. Figure 51a 
corresponds to the case where E > 0.0; Fig. 51b to E = 0.0; and Fig. 51c 
to E < 0.0. Clearly, the phase angle between the lift force and the 
velocity is smaller for the case of E > 0.0 than for E = 0.0. In fact, 

one can state that the said phase angle increases from a plus to a minus 

value. Only in the case of Fig. 51b is the phase angle such that the 

energy imbalance is for all intents and purposes zero. Furthermore, the 

stability of the system to small perturbations in amplitude is such that 
a small increase in amplitude leads to a decrease in phase between lift 
and displacement and thus to an increase in phase between lift and 
velocity. This, as noted earlier, leads to a decrease in amplitude. The 
reverse is true for a slight decrease in amplitude from a steady state 
osci I I at ion. 

3. Parametric Study of Synchronization 

Numerical experiments have been performed for the C and a^ combi- 
nations given in Table II in addition to that already presented. The 
results for each set are presented in Figs. 52 through 54. All of the 
parameters shown in these figures display similar variations with or 
U^. Furthermore, they are quite similar, at least in general form, to 
those experimentally obtained by Feng (see Fig. 47). A closer look at 
the experimental and calculated values show some differences, part of 
which may be attributable to the limitations of the numerical model and 
part to the limitations of the experimentation. In other words, all of 
the differences to be pointed out are not entirely consequences of the 
numerical model. 



141 




Fig. 51a. Variation of n, n» il and C|^ versus time for E>0.0. 




Fig. 51b. Variation of ri, ri,fi and C|^ versus time for E=0.0. 
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Fig. 51c. Variation of n> and C|^ versus time for E<0.0. 
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A comparison of Figs. 47 and 52 shows that the numerical model 

predicts a peak lift amplitude at a value smaller than U*. In Feng’s 

experiments, the lift and oscillation amplitudes reached their peak at 

nearly Identical U values. The reasons for this are not entirely clear. 

r 

It may be conjectured that it is a consequence of the extremely low damp- 
ing encountered in Feng's experiments = 0.001). In fact, mechanical 
systems with very small damping are much more finely tuned and more pre- 
cariously balanced than those with larger internal damping. It is there- 
fore natural to expect that the more stable systems experimented with in 
the present numerical calculations should exhibit larger differences in 
values at which they reach their extremals. 

The ratio of the cylinder oscillation frequency to the natural 
0) 0) 



C V 

frequency, — - — , remains practically constant at a value slightly under 
‘^n ^n 

unity in the experiments reported by Feng and others. In the numerical 

results, 0 ) /o) gradually increases from about 0.95 to 0.99 and then 
c n 

rapidly increases toward o) /w . The differences between the measurements 

o 

and calculations may be due to the difficulty of modelling some secondary 
features of flow which adjust themselves to render practically 

constant. 

The amplitudes of the lift coefficient and cylinder response cannot 
be compared with the experimental data because of the large differences 
in damping and mass parameters. However, it is interesting to note that 
the ratio of the maximum lift coefficient for the oscillating cylinders 
to the lift coefficient for the stationary cylinder (standard run) varies 
from 1.3 to 1.6 (see Figs. 48 and 52-54). This lift amplification is in 
conformity with most of the data reported in the literature [ 15 }, The 
cylinder response shows a dramatic drop for a small change in o)^ or U^, 
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remfniscent of the hysteritic drop in Feng's data. It should, however, be 
emphasized that the drop in n In the numerical calculations has nothing to 
do with the hysteresis resulting from the change in damping of the cylin- 
der. The sudden drop is a consequence of the changes in behavior of 
flow external to the cylinder. One further distinction must be made between 
the experimentally obtained sudden drop and that obtained numerically. In 
the experiments, the velocity of the flow (or U^) is increased gradually 
through the value at which the sudden change occurs. In the numerical 
calculations, or o)^ is kept constant at each point. Consequently, what 
is observed is not a sudden change in time but a change in amplitude as 

a function of o) . 

o 

The cylinder response for o)^ > w* does not remain at a constant 
amplitude. It is because of this fact that part of the response curves 
corresponding to that region are shown by broken lines in Figs. 48 and 
52-54. A sample plot of r) for > w* is shown in Fig. 55. 

In discussing the material damping and mass parameters, another 
parameter denoted by Sg = has been introduced. Among others, Vickery 

and Watkins C803 have attempted to show that q* (rij^ at perfect synchroni- 
zation) is uniquely determined by the value of a response or stability 
parameter, Sg. Sarpkaya C8l3 has shown that q* may be' fairly accurately 
represented by 

C, 

L o 

njl, = - (54) 

/0.06 + Sg2 

where C|^ is the lift coefficient for the oscillating cylinder. The values 
of C and a were judiciously selected to arrive at identical values of 
for three sets of calculations and twice the Sp value for one set of 
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Fig. 55. Representative variation of cylinder displacement versus time for 



calculations. According to the hypothesis advanced by Vickery and Watkins 
should be uniquely determinable by Sg. The results obtained through 
use of the model are tabulated below. The results in Table ill show that 
cylinder response at perfect synchronization is not very sensitive to the 
individual variations of the damping and mass ratios. 

TABLE I I I 



SET 




a 

o 




"M 


From Eq. (54) 


Numerical Results 


Eq. (56) 


I 


0.02 


0.0125 


1.6 


0.52 


. 0.54 


0.56 


1 1 


0.04 


0.0250 


1.6 


0.58 


0.58 


0.58 


1 1 1 


0.08 


0.050 


1 .6 


0.52 


0.61 


0.58 


IV 


0.04 


0.0125 


3.2 


0.24 


0.25 


0.26 



An exact relationship between Sg and other parameters may be 
established through the use of the energy balance at perfect synchroniza- 
tion. Re-writing Eq. (52) 

ttT t+T 




27t^S 



o 





Cj^ndt 



(52 repeat) 



and expressing ^ and as 

n = Dm sin w*t , C, = C sin (o)*t + 

M L L 

where 4>* is the phase angle between the force and the displacement at 
perfect synchronization one has 
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t+T 



J Cj^ndt 



= -^ = TTW *S 



o o t+T 



/ ^ 



Trails C, *sind)* 
o o L 

03*nJ 



n^dt 



(55) 



or 




Treats C, *s i n<f>* 
o o L ^ 

— — 



(56) 



Equation (56) shows that ri* is not a unique function of Sg. However, 
since do not vary significantly, varies approxi- 

mately as Sg * . A comparison of Eqns. (54) and (56) shows that Eq. (56) 
requires the use of parameters which must be determined either numerically 

or experimentally whereas Eq. (54) requires only a knowledge of and 

o 

0 )^ or U . 
o r 

The results obtained through the use of Eq. (56) for the four 
sets of input parameters are also shown in Table III. Clearly the direct 
and indirect predictions of n* are in excellent agreement, as would be 
expected. 

The n^-Sg relationship given by Eq. (56) is expected to hold true 

for values of o) other than w* provided oa is confined to the region of 
o o o 

synchronization. Any small deviations between the predictions of 



TToj S C, s i nd) 
_ _ o o L 

" 0) 



(56a) 



from those directly calculated with the model are simply due to small 
modulations in the amplitude of oscillations. Numerical calculations which 
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can easily be performed through the use of Figs. 48 and 52-54 have shown 
that the predictions of Eq.(56a) are within a few percent of those obtained 
directly. 

The transverse oscillations of a cylinder in uniform flow Increase 
the affective projected area of the cylinder and hence the bluff ness of 
the body. Consequently, hydroelastic oscillations of a cylinder, such as 
the strumming of a cable, not only cause alternating stress, fatigue, 
and noise but also increased drag. 

The numerical model calculated two important characteristics of 

the fluctuating drag force: the mean drag and the amplitude of oscillation. 

Figure 56 shows the mean drag coefficient, Cp, and the amplitude of the 

fluctuating component, ACp, as a function of or U^. Also shown in 

Fig. 56 i s a data point obtained experimentally CSlJ.for the corresponding 

( 1 ) and q,, values, 
o M 

Figure 57 i s a representative plot of the instantaneous value of 
Cp as a function of time. It shows, among other things, that the drag 
force oscillates at a frequency twice that of the lift force (for lift 
T 10 and for drag T = 5). 

The separation points on an oscillating cylinder oscillate with 
a frequency equal to that of the cylinder. Representative results pre- 
sented in Fig. 58 show that the mean position of the separation angle 
remains practically constant at a value of 0^ = 77 degrees. The amplitude 
of the fluctuations follow a pattern very similar to that of rij^. The 
stagnation point undergoes similar oscillations. The results have shown 
that the instantaneous position of the stagnation point nearly coincides 
with that which would have' resulted with a uniform flow at the instantaneous 
total velocity. In other words, the small oscillations of the stagnation 
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Fig. 56. Drag coefficient versus for C=0.02, a^=O.OI25. 
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Fig. 58. Representative variation of separation angle versus time for synchronized 
oscillations (^=0.02, a =0.0125, (o =1.125). 



point relative to the uniform flow are obscured by the larger oscillations 
resulting from the oscillations of the cylinder or of the ambient flow 
relative to the cylinder. 
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VII. CONCLUSIONS 



The results presented herein warrant the following conclusions: 

1. A discrete-vortex model based on the rediscretization of shear layers 
and the use of appropriate vorticity cancellation mechanisms can produce 
results which are in conformity with those obtained experimentally; 

2. The numerical model allows experimentation with the variation of any 
one of the independent parameters and permits one to understand the 
governing physical mechanisms; 

3. It has been shown that two major flow parameters, namely the drag co- 
efficient and the Strouhal number, are least sensitive to changes in the 
characteristics of the wake. The Strouhal number is intrinsic to the 
shape of the bluff body and its magnitude cannot be varied more than a 
few percent by varying the strength of the vortices by physical means. 
Experiments over the past 50 years are in conformity with these conclu- 
sions; 

4. The lift coefficient is the one most sensitive to the changes of 
flow in the near wake. This finding of the numerical model and the 
difficulty of obtaining consistent laboratory data for lift are in con- 
formity with each other. The numerical model explains the physical 
mechanism as to why lift should be a strong function of the evolution 
of the near wake; 

5. The discrete vortex model, as all other numerical methods attempting 
to model a turbulent region of a flow and its interaction with a laminar 
region, is in need of a hypothesis to properly simulate the circulation 
dissipation and the vorticity diffusion caused by turbulence. 
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Notwfthstandfng this fact, a very simple decay mechanism is capable of 
simulating the said consequences of turbulence; 

6. The predictive powers of the numerical model have been put to test 
through its use in the study of self-excited oscillations of elastically 
mounted circular cylinders. The results were found to be in conformity 
with practically all observations and measurements made in recent years 
and led to the understanding of the mechanisms governing the synchroni- 
zation phenomenon. 
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